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ABSTRACT 

Aims. Our goal is to identify bipolar H ii regions and to understand their morphology, their evolution, and the role they play in the 
formation of new generations of stars. 

Methods. We use the Spitzer-GLIMPSE, -MIPSGAL, and Herschel-Ki-GAh surveys to identify bipolar Hii regions, looking for 
(ionized) lobes extending perpendicular to dense filamentary structures. We search for their exciting star(s) and estimate their distances 
using near-IR data from the 2MASS or UKIDSS surveys. Dense molecular clumps are detected using Herschel-SFIRE data, and we 
estimate their temperature, column density, mass, and density. MALT90 observations allow us to ascertain their association with the 
central H ii region (association based on similar velocities). We identify Class O/I young stellar objects (YSOs) using their Spitzer and 
Herschel-FACS emissions. These methods will be applied to the entire sample of candidate bipolar Hii regions to be presented in a 
forthcoming paper. 

Results. This paper focuses on two bipolar Hii regions, one that is especially interesting in terms of its morphology, G319.88+00.79, 
and one in terms of its star formation, GOlO.32-00.15. Their exciting clusters are identified and their photometric distances estimated 
to be 2.6 kpc and 1.75 kpc, respectively; thus GOlO.32-00.15 (known as W31 north) lies much closer than previously assumed. We 
suggest that these regions formed in dense and flat structures that contain filaments. They have a central ionized region and ionized 
lobes extending perpendicular to the parental cloud. The remains of the parental cloud appear as dense (more than 10"^ cm"^) and 
cold (14-17 K) condensations. The dust in the photodissociation regions (in regions adjacent to the ionized gas) is warm (19-25 K). 
Dense massive clumps are present around the central ionized region. GOlO.32-00.14 is especially remarkable because five clumps of 
several hundred solar masses surround the central Hii region; their peak column density is a few 10^^ cm“^, and the mean density in 
their central regions reaches several 10^ cm"^. Four of them contain at least one massive YSO (including an ultracompact H ii region 
and a high-luminosity Class I YSO); these clumps also contain extended green objects (EGOs) and Class II methanol masers. This 
morphology suggests that the formation of a second generation of massive stars has been triggered by the central bipolar H ii region. 
It occurs in the compressed material of the parental cloud. 

Key words. ISM: Hii regions - ISM: dust - ISM: individual objects: GOlO.32-00.15, G319.88+00.79 - Stars: formation 


1. Introduction 

Hi-GAL, the Herschel infrared GALactic Plane survey (Molinari 
et al. 2010a), gives an overall view of the distribution of the 
different phases of the Galactic interstellar medium (ISM). 
According to Molinari et al. (2010b), “The outstanding feature 
emerging from the first images is the impressive and ubiqui¬ 
tous ISM filamentary nature.” Now that all of the Galactic plane 
has been surveyed by Hi-GAL we still have the same view of a 
very filamentary structure, especially at HerschelSVlE^ wave¬ 
lengths. Several papers present a preliminary analysis of fila¬ 
ments (for example, Arzoumanian et al. 2011, 2013; Hill et al. 
2011; Peretto et al. 2013) and suggest that the filaments are the 
main birth site of prestellar cores (for example, Molinari et al. 
2010b; Andre et al. 2010). 
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The Spitzer-GLlMPSE survey has shown a “bubbling 
Galactic plane” (Churchwell et al. 2006, 2007; Simpson et al. 
2012). These bubbles observed at 8.0 pm enclose classical Hii 
regions (Deharveng et al. 2010; Anderson et al. 2011). The ex¬ 
act morphology of most of these bubbles is still unknown; are 
they 3D spherical structures or 2D rings? The vicinity of 43 
bubbles has been observed by Beaumont & Williams (2010) in 
the (3-2) transition. These observations suggest that these 
bubbles are not spherical structures, but rings, and that Hii re¬ 
gions form in flat molecular clouds (with a thickness of a few 
parsecs). If massive stars form in filamentary ID or 2D struc¬ 
tures, thus dense filaments or sheets, so in a medium where den¬ 
sity gradients are present, many bipolar nebulae should form and 
be observed. The Green Bank Telescope Hii Region Discovery 
Survey (Bania et al. 2010; resolution 82") shows the opposite 
situation, the rarity of bipolar nebulae, hence the conclusion of 
Anderson et al. (2011), based on the comparison of the radio 
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and Spitzer 8.0 pm emissions and on the radio recombination 
linewidths, that the majority of bubble sources are 3D structures. 
Statistics on bipolar nebulae should resolve this controversy. 

Bipolar Hii regions are also interesting because they have 
a simple morphology, allowing one to locate the ionized and 
neutral components in space in their associated complex. 
And, as predicted by Fukuda & Hanawa (2000; Sect. 3), star 
formation can be triggered by the expansion of such nebulae 
adjacent to dense filaments. In the following, we use Spitzer and 
Herschel observations of the Galactic plane to search for these 
nebulae. 

Using the Spitzer GLIMPSE survey, and the Herschel Hi- 
GAL survey we have identified 16 candidate bipolar Hii regions 
in a zone of the Galactic plane between +60° in longitude and 
±1° in latitude. In this first paper we describe two of them, 
G319.88+00.79 and GOlO.32-00.15, that we consider the most 
outstanding in terms of morphology and star formation, respec¬ 
tively. In this paper, we give in Section 2 a quick description 
of the observations we use for this study, mainly the Spitzer- 
GLIMPSE and -MIPSGAL surveys, the Herschel Hi-GAL 
survey, and the MALT90 survey. The few existing models of 
bipolar nebulae are discussed in Section 3, and we discuss how 
these nebulae can be identified and how star formation in their 
vicinity can be observed. The general methods used to estimate 
the physical parameters of various structures (molecular clumps, 
cores, etc.) are described in Section 4, where we also explain 
how we detect young stellar objects (YSOs), especially candi¬ 
dates Class O/I YSOs. Sections 5 and 6 are devoted to the study 
of the G319.88+00.79 and GOlO.32-00.15 bipolar nebulae. We 
discuss our findings in Sect. 7 and give conclusions in Sect. 8. 
A second paper will contain the catalogue of bipolar nebulae, 
where a study of each of them uses the same observations and 
the same methods as described in this first paper, and our gen¬ 
eral conclusions concerning star formation in their surroundings. 

Eollowing Eoster et al. (2011) we refer to molecular 
“clumps” when speaking about molecular structures with typical 
masses in the range 50 M© - 500 M©, typical sizes in the range 
0.3 pc - 3 pc, and typical mean densities in the range 10^ cm“^ 
- 10^ cm“^. These structures may contain smaller (<0.1 pc) and 
denser (>10"^ cm"^) substructures called “cores”. 

In the following we use the term “filament” to describe a 
cylindrical structure (ID), and “sheet” to describe a flat struc¬ 
ture (2D). It has a small thickness in one direction and a larger 
section perpendicular to this direction. A sheet observed edge- 
on appears as a filament (and so some structures called filaments 
may be sheets). 

2, Observations 

We use several surveys for this study. We summarize their 
wavelengths and resolutions in Table 1. The Hi-GAL survey 
(Molinari et al. 2010a) covers a two-degree-wide strip of the 
Galactic plane, in five bands centred at 70 pm, 160 pm, 250 pm, 
350 pm, and 500 pm. The data reduction pipeline is described 
by Traficante et al (2011). The Hi-GAL maps, which trace the 
dust continuum emission, are complemented by: 

• DSS2-red or SuperCOSMOS Ha maps (Parker et al. 2005) 
and radio-continuum maps (NVSS survey, Condon et al. 1998; 
SUMSS survey, Bock et al. 1999) to describe the emission of the 
ionized gas (respectively Ha emisssion or free-free emission). 
We use in particular the MAGPIS survey (Multi-Array Galactic 
Plane Imaging Survey; Helfand et al. 2006) and the CORNISH 


survey (Hoare et al. 2012; Purcell et al. 2013) to detect ultracom¬ 
pact (UC) Hii regions. Radio recombination lines are also used 
to give the velocity of the ionized gas. 

• 2MASS (Skrutskie et al. 2006), UKIDSS (from the 
Galactic plane Survey; Lawrence et al. 2007), or VISTA (Minniti 
et al. 2010) maps and catalogues to describe the stellar content 
of the selected regions. 

• Spitzer-GIAMPSE maps at 3.6 pm, 4.5 pm, 5.8 pm, 8.0 pm 
(Benjamin et al. 2003), and Spitzer-MlFSGAh maps at 24 pm 
(Carey et al. 2009) to detect candidate YSOs, and describe the 
distribution of specific dust grains. WISE (Wide-Eield Infrared 
Survey Explorer; Wright et al. 2010) maps are used if the Spitzer 
maps are missing or saturated. The Spitzer-GLlMPSE and - 
MIPSGAL catalogues are used, as well as the Robitaille et al. 
(2008) catalogue of intrinsically red sources. 

• Maser observations of various molecular species (espe¬ 
cially methanol) and observations of extended green objects 
(EGOs; Cyganowski et al. 2008) to indicate the presence of 
YSOs with outflows (Breen et al. 2010, and references therein). 

• MALT90 (Millimeter Astronomy Legacy Team 90 GHz) 
observations of the molecular clumps present in the vicinity of 
the Hii regions. A general description of the MALT90 survey is 
given by Jackson et al. (2013). About 2000 dust clumps observed 
at 870 pm (ATLASGAL survey, Schuller et al. 2009) have been 
mapped with the Mopra telescope at a frequency ~90 GHz and a 
beam size of 38". Each map has a size of 3'x3', with a pixel size 
of 9". Sixteen molecular lines are mapped simultaneously with 
a spectral resolution of 0.11 km s“^ More details can be found 
in Eoster et al. (2011, 2013). Here we use these observations to 
obtain the velocity fields of the molecular clumps, and ascertain 
their association with the ionized regions. 


Table 1. List of the main surveys used for this study 


IR surveys 

Wavelength 

FWHM 


Oum) 

(") 

2MASS J,H,K, 

1.25, 1.65,2.16 

>2.5 

UKIDSS J,H,K 

1.25, 1.63,2.20 

<1 

\ISTAJ,H,K, 

1.25, 1.65,2.15 

<1 

Spitzer-IRAC 

3.6, 4.5, 5.8, 8.0 

<2 

WISE 

3.4, 4.6, 11,22 

6.1, 6.4, 6.5, 12.0 

Spitzer-MIFS 

24 

6 

Hi-GAL 

70, 160, 250, 350, 500 

8.5, 12.5, 18,25,36 

Radio surveys 

Wavelength 

FWHM 


(cm) 

(") 

NVSS 

21 

45 

SUMSS 

35 

45x45eosec(5 

MAGPIS 

20 

6 

CORNISH 

6 

1.5 


3. Simulations of bipolar nebulae and triggered star 
formation 

Very few studies describe the formation and evolution of a 
bipolar Hii region^. Bodenheimer et al. (1979) present a 2D 
hydrodynamic simulation following the evolution of an Hii 
region excited by a star lying in the symmetry plane of a flat 
homogeneous molecular cloud surrounded by a low-density 

Tn contrast, many studies consider the formation of bipolar plane¬ 
tary nebulae, but the situation is very different sinee these result from 
the ejection of a small amount of material by an evolved low-mass star. 
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medium. (The parental cloud is in pressure equilibrium with 
the surrounding medium.) Figure 1 illustrates the formation 
of a bipolar Hii region according to this simulation. Initially 
(schema a), a spherical Hii region grows inside the cloud. A 
bipolar nebula forms when the ionization front breaks through 
the two opposite faces of the cloud simultaneously (schema 
b). Then, the dense ionized material pours out of the dense 
cloud with velocities of up to 30 km s“\ forming a double cone 
structure (schema c); with time the angle of the cone widens, 
and the density decreases in the central parts of the Hii region; 
neutral material accumulates at the waist of the bipolar nebula. 
These H ii regions appear bipolar when the angle between the 
line of sight and the cloud’s symmetry plane is small (when the 
flat parental cloud is almost seen edge-on). 

Fukuda & Hanawa (2000) study sequential star formation 
in a filament, triggered by the expansion of a nearby Hii re¬ 
gion. Their 3D simulations take both (magneto) dynamical and 
gravitational effects that may induce star formation into account. 
In these simulations (Fig. 2) a dense filamentary cloud is com¬ 
pressed by a nearby expanding H ii region and fragments to form 
cores sequentially. The H ii region is initially off the filament; as 
it expands and grows in size it interacts with the filament. The 
filament is compressed and pinched, and it possibly separates 
into two parts^ ; the H ii region appears bipolar to an observer if 
the line of sight is roughly perpendicular to the plane contain¬ 
ing the filament and the exciting star. The dynamical compres¬ 
sion of the filament triggers the formation of a first-generation 
of cores at the waist of the bipolar nebula. With time the density 
in the cores increases (from compression by the H ii region and 
self-gravity), they become gravitationally bound and they col¬ 
lapse. Star formation probably occurs, but this is not followed 
in the Fukuda & Hanawa paper. Later, the formation of the first- 
generation cores changes the nearby gravitational held and in¬ 
duces the formation of second-generation cores. The time of the 
formation of the first-generation cores depends upon the distance 
of the exciting star to the axis of the filament, of the energy of the 
expansion (through its velocity and the radius of the H ii region), 
and of the magnetic held. Depending on the distance between 
the exciting star of the Hii region and the filament’s axis, the 
number of first-generation cores ranges from 2 to 4. 

Fukuda & Hanawa did not simulate what happens if an H ii 
region forms in or near a dense sheet of material (in a flat struc¬ 
ture instead of a filament). 

3.1. What do we expect to see? 

The signatures of bipolar Hii regions that we expect to see are: 

• A dense filament, observed as a cold elongated high col¬ 
umn density feature. This structure will be observed in emission 
at Herschel-SFIRE wavelengths. It can appear in absorption, as 
an elongated infrared dark cloud (IRDC) at 8.0 pm or 24 pm. 
This “filament” could also be a sheet seen edge-on. 

• An HII region, bright in its central part and displaying two 
ionized lobes perpendicular to the filament. It will appear as 
a radio-continuum source, centred on the filament and slightly 
elongated in two directions perpendicular to it. The two lobes 
will be especially well traced at 8.0 pm, a wavelength band dom¬ 
inated by the emission of PAHs located in the photodissociation 


^Fukuda & Hanawa (2000) only represent in their figures the dis¬ 
tribution and kinematics of the neutral material. No limits are given for 
the ionized region. The limits drawn in Fig. 2 correspond to our best 
guess. 


region (PDR) surrounding the ionized gas and excited by the 
UV radiation leaking from the H ii region (PAHs are destroyed 
inside the ionized region, Povich et al. 2007; Pavlyuchenkov et 
al. 2013). At 24 pm we will see an extended central source, 
similar to the radio-continuum source (Deharveng et al. 2010; 
Anderson et al. 2011); also 24 pm emission will be observed in 
the PDRs surrounding the two ionized lobes. The 24 pm emis¬ 
sion observed in the direction of the ionized region comes from 
very small dust grains (size of a few nanometres) located inside 
the ionized gas and out of thermal equilibrium after absorption of 
ionizing photons (Pavlyuchenkov et al. 2013). At 70 pm a small 
fraction of the emission comes from the ionized region, but the 
bulk of the emission comes from the PDRs surrounding the ion¬ 
ized region (emission of silicate grains; size of a few hundred 
nanometres; Pavlyuchenkov et al. 2013). 

• The exciting star(s) or cluster. If the extinction is not too 
high in its direction, it should be observed in the near IR, near 
the symmetry axis of the filament, at the centre of the H ii region. 

• Dense material and possibly dense cores at the waist of the 
bipolar nebula, containing YSOs or not, depending on whether 
star formation is presently at work or not. If present. Class O/I 
YSOs will possibly appear as point-like sources at Herschel 
wavelengths, especially at 70 pm'. Class I/II YSOs will appear 
as Spitzer-mi& IR sources, prestellar cores will be without any 
associated mid-IR sources. 

• More cores located along the filament, but not adjacent to 
the H II region, if second-generation cores have already formed 
(according to the simulations of Fukuda & Hanawa 2000). If 
present, they will appear as HerschelS?\RE compact sources. 

Of course, all the components of the complex, the filament, 
the clumps, and the Hii region must have similar velocities to 
ascertain their association. 

In the simulation of Bodenheimer et al. (1979) the medium 
surrounding the molecular condensation has a very low density. 
Thus, it is very quickly ionized, and the two lobes are density 
bounded and not bounded by ionization. If this medium is denser 
(but however less dense than the parental cloud), we can expect 
the two lobes to be ionization-bounded and surrounded by an 
ionization front (IF) of D type^. Then they should appear as 
closed lobes surrounded by PAHs emission, and also possibly 
surrounded by neutral material collected during the expansion 
of the ionized gas. Thus the presence of a layer of neutral mate¬ 
rial surrounding the ionized lobes can also be a feature of bipolar 
nebulae. 

We search for these signatures and discuss them for each 
candidate bipolar H ii region. 

Several objects can be mis-identifled as bipolar Hii regions. 
They are: 

• Nebulae associated with evolved stars, such as G308.7- 
00.5 (Fig. 3). These “well defined nebulae are formed during a 
brief post-main-sequence phase when a massive star becomes a 
Wolf-Ray et star or a Luminous Blue Variable star” (Gvaramadze 
et al. 2010). Some of them show a bipolar structure when ob¬ 
served by Spitzer. 

• Adjacent bubbles surrounding distinct H ii regions, such as 
NIO and Nil (Fig. 3). According to the Churchwell et al. cat¬ 
alogue (2006) “BP indicates a bipolar bubble or a double bub¬ 
ble whose lobes are in contact.” NIO and Nil fall within the 


^An IF is of D type when the expansion velocity of the Hii region 
is comparable to the sound speed in the ionized gas, of the order of 
10 km s ~^. This IF is preceded on the neutral side by a shock front; with 
time neutral material accumulates between the two fronts. 


3 



L. Deharveng et al.: Bipolar nebulae 



Fig. 1. Formation of a bipolar H ii region, according to the simulation of Bodenheimer et al. (1979; their figure 4; the thickness of the parental plane 
is 1.3 pc, its density 300 H 2 cm"^). The ionized material appears in red, the neutral material in blue; the density of the material is indicated by the 
saturation level. Schema a shows the dense molecular plane surrounded by low density material, the initial Stromgren sphere around its central 
exciting star (T^=4xl0'^ K). In schema b, at 3x10"^ yrs, the expanding Hii region reaches the border of the parental cloud, and the low density 
material is quickly ionized; a bipolar nebula forms. In schema c, at IxIO^ yrs, the high density ionized material flows away from the central region; 
high density molecular material accumulates at the waist of the bipolar nebula, forming a torus of compressed material. 



Fig. 2. Star formation triggered by the expansion of an H ii region close to a filament, according to the simulation of Fukuda & Hanawa (2000; 
their figure 1; model Cl/’: neutral sound speed 0.3 km s“^ maximum density on the axis of the filament 2x10^ cm“^, exciting star at 0.025 pc 
of the filament’s axis). The neutral filament appears in grey, the ionized gas in pink; we show what we believe to be the extent of the ionized gas 
in the plane containing the axis of the filament and the exciting star (black cross) of the H ii region. Schema a shows the expanding H ii region 
(age 8.5x10"^ yr) pinching the molecular filament. In schema b, first-generation cores have formed by compression of the neutral material by the 
ionized gas (age 2.1x10^ yr). In schema c, second-generation cores have formed due to gravity (age 3.0x10^ yr). The first- and second-generation 
cores are separated by O.I pc. 


last case. Figure 3 shows that the 8.0 jim bubbles enclose two 
separate Hii regions: we clearly see two separate zones of ex¬ 
tended 24 yum emission at the centre of each bubbles. In bipolar 
Hii regions like Sh 201 (see below), the extended 24 yum emis¬ 
sion comes from the central region, the region where the lobes 
originate. Furthermore, the HerschelSV^RE maps show no cold 
filament between NIO and Nil. 

The Hii region Sh 201 is a textbook example of a bipolar 
Hii region (Fig. 3). It has been studied using deep Herschel 
observations (see Appendix C of Deharveng et al. 2012). The 
complex displays a cold filament observed at SPIRE wave¬ 
lengths, an HII region centred on this filament, and two lobes, 
well traced at 8.0 yum (Fig. 3), perpendicular to the filament. 
The ionized region is bright in its central part. It shows a narrow 
waist, limited on each side by two massive clumps. YSOs of 
Class 0 and Class I are observed in the direction of these clumps 


(Deharveng et al. 2012). 

In all the figures of this paper the regions are presented in 
Galactic coordinates. 


4. Reduction methods 

4.1. Determination of the physical parameters of the 
filaments and clumps 

For all estimates in this paper, we assume a gas-to-dust ratio 
of 100. We assume the following values for the dust opacity 
Ky, 35.2, 14.4, 7.3, and 3.6 cm^ g“\ respectively, at 160, 250, 
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Fig. 3. Composite colour image of H ii regions, as seen by Spitzer. Red is the MIPSGAL image at 24 pm showing the emission of the hot dust, 
green is the GLIMPSE image at 8.0 pm dominated by the PAHs emission from the PDRs, and blue is the GLIMPSE image at 4.5 pm showing 
the stellar sources. Sh 201 is a textbook example of a bipolar Hii region; the red column density contours trace the parental molecular filament 
(levels of 2, 5, and 15 xlO^^ cm"^). NIO and Nil are adjacent bubbles around distinct Hii regions. This is a case of mis-identified bipolar nebula. 
G308.7-00.5 is a nebula associated with an evolved star. 


350, and 500 The opacities are uncertain, as is the 

gas-to-dust ratio. The dust opacity probably differs depending 
on the content of the clumps, for example whether they have 
internal sources (Paradis et al. 2014). This remains the main 
source of uncertainty of the temperature, column density, and 
mass determinations. (The column density and the mass are 


"^These values correspond to the formula 

which assumes a spectral index j3 of 2 for the dust emissivity and gives 
a dust opacity at 300 pm of 10.0 cm^ g"^ (Andre et al. 2010). These 
values have been discussed by Sadavoy et al. (2013); they show that as¬ 
suming [5-2 with the Herschel bands alone (using the 160 pm - 500 pm 
range) affects the measured best-fit temperature by <2 K if p- 1.5 or 
= 2.5, and that “the Herschel-only bands provide a decent first look at 
the clump masses for an assumed value ofp- 2”. 

^The dust opacities determined by Ossenkopf & Henning (1994; 
dust with thin ice mantle, density 10^ cm“^) are also often used; they 
correspond to a ~1.8 and have values of 40, 18, 10, and 5 cm^ g"^ 
respectively at 160 pm, 250 pm, 350 pm, and 500 pm. 


inversely proportional to /Cy.) 

We have obtained dust temperature maps for all the regions, 
using Herschel data between 160 pm and 500 pm. We assume 
that the dust emission is optically thin at these wavelengths. We 
do not use the 70 pm data; the 70 pm maps show that a fraction 
of this emission comes from ionized regions, as is the case for a 
larger portion of the 24 pm emission; i.e. the 24 pm and 70 pm 
emissions probe a different gas/dust component than emission 
observed at longer wavelengths (Sect. 3.1). Furthermore, some 
of the clumps are not optically thin at 70 pm', for example, most 
of the clumps in the field of GOIO.32-00.15 (Table 4). The tem¬ 
peratures that we determine are those of cold dust grains in ther¬ 
mal equilibrium, located in the PDRs of H ii regions and in the 
surrounding medium. The steps in the temperature determina¬ 
tion are: 1) the 160 pm, 250 pm, and 350 pm maps have been 
smoothed to the resolution of the 500 pm map; 2) all the maps 
have been regridded so they have the same pixel size of 11.5" 
as the 500 pm map (at the same location); 3) the SED of each 
pixel has been fitted, using a modified blackbody model; thus a 
first temperature map. We applied colour corrections to the tem- 
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perature map of the G319.88+00.79 field (using the first tem¬ 
perature map, and following the PACS and SPIRE documenta¬ 
tion). The colour-corrected temperatures are higher than those 
obtained without corrections, with a mean difference of 0.58 K 
(in the range 0.3 K - 0.7 K). However, since the colour correc¬ 
tions were still uncertain at the time we created the temperature 
maps, and because they are not very large (less than 6% on the 
fluxes,) we decided to present in the following the temperature 
maps obtained without colour-corrections. 

No background subtraction has been made, so the tempera¬ 
ture obtained is a mean temperature for the dust along the line of 
sight, weighted by the emission. If assigned to a specific struc¬ 
ture, it assumes that the emission from the structure is dominant, 
which is possibly not the case. These temperature maps will be 
used to study the temperatures of different structures present in 
the field. 


The molecular hydrogen colum density is given by 

F 10^^ 

A^(H2) = 100 J' — 

Ky 5y(T(5yst) 2.8 


( 2 ) 


where the brightness Fy is expressed in MJy sr“^ (the units of the 
Hi-GAL maps), the Planck function By in J s“^ m“^ Hz“\ the 
factor of 2.8 is the mean molecular weight, mn is the hydrogen 
atom mass. 

To keep the best possible spatial resolution, we use temper¬ 
ature and column density maps regridded the resolution of the 
250 yum data. The original temperature maps have the resolution 
of the 500 Jim maps; however, since the dust temperature varies 
only smoothly across the fields, we expect that this extrapolated 
temperature does not differ too strongly from the one that would 
be obtained if all the maps had the resolution of the 250 jim 
observations. Here again, no background emission has been 
subtracted, and thus, the column density corresponds to the 
whole line of sight. 


Two methods can be used to estimate the mass of a structure. 
In the first method, we integrate the column density map within 
an aperture. To compare the parameters of the clumps observed 
in the vicinity of the bipolar nebulae, we use apertures defined 
by the intensity at half the peak value. We prefered to use this 
column density method for extended structures. In the second 
method, after assuming a mean temperature for the structure (ei¬ 
ther by fitting the SED of the structure or by integrating the tem¬ 
perature map), we determine the mass of the structure following 
the method of Hildebrand (1983). The total (gas+dust) mass of 
a feature is related to its integrated flux density Sy by 


^(gas+dust) — 


Sy 

Ky Sy(7^dust) 


( 3 ) 


where D is the distance to the source and 5y(rdust) is the Planck 
function for a temperature T^ust- Here again we use the 250 jim 
fluxes (to have the best resolution). This method is preferably 
used for point sources. We have used both methods to esti¬ 
mate the mass of the three clumps in the G319.88+00.79 field 
(Table 2). The masses obtained differ by less than 20% (see 
Table 2). 

If a clump can be modelled as an elliptical Gaussian of uni¬ 
form temperature, its total mass is twice the mass measured us¬ 
ing an aperture following the level at half the peak’s value. Also, 
the derived mass is not dependent on the beam size, which there¬ 
fore allows us to compare the masses of clumps at different dis¬ 
tances and angular sizes. 


We also estimate the mean density in the central regions of 
the clumps, using their mass and beam deconvolved size ob¬ 
tained with an aperture following the level at half the maximum 
value, assuming a spherical morphology and a homogeneous 
medium. 

A background subtraction is done in only one instance: to es¬ 
timate the parameters of the bright clumps discussed in Sects. 5.3 
and 6.3 more accurately (details are given in the text and in the 
notes to Tables 2 and 4). 


4.2. The velocity field of the molecular clumps 

The MALT90 observations are used to obtain the velocity field 
of each molecular clump. In both regions we used the bright¬ 
est molecular transitions, H^^CO^ (1-0), N 2 H^ (1-0), and HNC 
(1-0). We performed Gaussian fittings of the main lines to re¬ 
trieve the distribution and velocity field of the molecular gas. 
We also used the H^^CO^ (1-0) to confirm the localized exis¬ 
tence of H^^CO^ self-absorption, and to validate the use of a 
single Gaussian fit of the double profile of H^^CO^ as providing 
a good approximation of the systemic velocities. 

We present the velocity fields resulting from this Gaussian 
fitting, and use those to ascertain the association of the observed 
molecular clumps with the HII regions. Details about the reduc¬ 
tion procedure can be found in Zavagno et al. (in preparation). 


4.3. Identification of young stellar objects 

In the field of each region we try to detect Class 0 or Class I 
YSOs and prestellar cores, if any. We do not consider Class II 
YSOs because, depending on the mass of their central source, 
they can possibly be as old as the exciting stars of the bipolar 
nebulae. Only Class O/I YSOs, with an age <10^ yr (Andre et 
al. 2000 and references therein) compared to an age of 10^ yr 
or more for extended Hii regions, are most probably second- 
generation forming stars. Different indicators have been used to 
identify candidate Class I or Class 0 YSOs: 

• The Spitzer [3.6]-[4.5] versus [5.8]-[8.0] diagram (Allen 
et al. 2004). We have used the Spitzer GLIMPSE Source 
Catalogue (which combines the three surveys GLIMPSE I, II, 
and 3D) to draw these diagrams, measuring a few missing mag¬ 
nitudes when necessary^. The limits in these diagrams of the 
location of Class I and Class II sources are only indicative, since 
the position of the sources depends of the external extinction, 
which is often unknown, and relies on a good background esti¬ 
mate, which is often difficult to obtain. 

• Emission of point sources at 24 pm. We used the 24 pm 
magnitudes from the Robitaille et al. (2008) catalogue of in¬ 
trinsically red objects anf from the MIPSGAL catalogue, when 
available; otherwise, we used these magnitudes to calibrate our 
own measurements. The spectral index, given by 


dlog(AF^) 

dlog(A) 


( 4 ) 


has been estimated between K and 24 pm, when possible, to 
confirm or estimate the nature of the point sources. According 
to Lada (1987) and Greene et al. (1994), Class I YSOs have 
a >0.3 in the range ^-20 pm; following Greene et al. sources 


^Some interesting sources seen in the Spitzer images are missing 
in the catalogues. Most often such sources are difficult to measure be¬ 
cause they are found against a bright background. Most such sources are 
observed in the direction of PDRs, and therefore are interesting when 
discussing star formation triggered by Hii regions. 
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with -0.3< a <0.3 are ”flat-spectmm” sources (of uncertain 
evolutionary status between Class I and Class II), and sources 
with a <-0.3 are Class II or Class III YSOs. This index is how¬ 
ever strongly affected by an external extinction^ that is generally 
unknown because we do not know the position along the line 
of sight of sources observed in the direction of clumps. Thus 
we prefer to consider the colour [8.0]-[24] of the point sources, 
since it is almost unaffected by extinction; a >0.3 corresponds 
to a colour [8.0]-[24]>3.9 mag, a <-0.3 to [8.0]-[24]<3.2. 

• Emission of compact sources at Herschel wavelengths - 
mainly the detection of 70 jim point sources, the counterparts 
of Class O/I YSOs or prestellar cores. Owing to the resolution 
of Herschel data and the distance of our regions, it is difficult 
to distinguish between point sources and compact sources in our 
fields. Most of our 70 yum point sources have no detectable coun¬ 
terpart at longer wavelengths. What can we deduce from their 
70 yum fluxes? 

Dunham et al. (2008) find a tight correlation between the 
70 yum flux and the internal luminosity of low-mass protostars 
(that of their central objects), and it is valid up to 50 L©, which 
corresponds to a flux of 16 Jy at 70 yum for a source at a distance 
of 1 kpc. Most of our measured sources are brighter at 70 yum 
(see Tables 3 and 5). Ragan et al. (2012) have shown (their fig¬ 
ure 8) that the 70 yum luminosity of more massive sources cor¬ 
relates well with their total luminosity. There is good agreement 
between the correlations established for low and high-luminosity 
sources. We use this correlation to roughly estimate the luminos¬ 
ity of sources detected only at 70 yum. 

Stutz et al. (2013) have shown that a ratio 
/^?g(TF^(70)/TF^(24))>1.65 defines what they call ‘TAGS 
Bright Red sources” (PBRs), which are extreme Class 0 sources 
with high infall rates. We estimate this ratio and identify such 
objects when possible. 

As discussed by Robitaille et al. (2008, and references 
therein), evolved stars (AGB stars) can be mistaken for YSOs. 
We have used the criteria recommended by Robitaille et 
al. (2008) to identify candidate xAGB stars (extreme AGB 
stars; AGBs with very high mass-loss rates): [4.5]<7.8 mag. 
When possible we measured the [8.0]-[24] colour, because 
another criterion to identify AGB stars is [8.0]-[24]<2.5 mag 
(Robitaille et al. 2008). However, as stressed by Robitaille et al., 
the separation between AGBs and YSOs is only approximate. 

The magnitudes of Spitzer and Herschel point sources are 
measured, if necessary, using aperture photometry for isolated 
sources superimposed on a faint background or using PSF fitting 
with DAOPHOT (Stetson 1987) in more difficult cases (sources 
superimposed on a bright variable background; details are given 
in Deharveng et al. 2012). 


5. G319.88+00.79 

This nebula is a textbook example of a bipolar H ii region. It is 
composed of the two Spitzer bubbles S97 and S98 (Churchwell 
et al. 2006). Figure A (a) (displaying the central 24 pm emission 
region) and Fig. 5 (c) (displaying the central radio-continuum 
emission region) show that the two bubbles only enclose one 
Hii region, bright in its central part. Figure A (h) shows that a 

^An uncorrected (visual) external extinction ~50 mag (correspond¬ 
ing to a column density ~5xl0^^ em"^) deereases the value of a by 
-1.04. 


bright cold dust filament is present at the waist of the nebula, 
which is conspicuous on all the Herschel SPIRE images (and 
even more conspicuous in the column density map. Fig. 7 (d)). 
The two lobes, perpendicular to the filament, are observed at 
8.0 yum, 24 yum, and 70 pm. All the region is symmetric with 
respect to the filament. However, one ionized lobe is enclosed 
well by the S98 bubble, whereas the other one (S97 bubble) is 
larger and open. Thus the H ii region is ionization-bounded on 
one side and possibly density-bounded on the other side. The 
SPIRE images shows that the lobes (especially the upper one) 
are surrounded by a low emission layer. This emission possibly 
comes from neutral material collected during the expansion of 
the ionized gas. 

G319.88+0.79 is an optical Hii region. Figures 5 and 6 show 
the central region at the waist of the bipolar nebula. We very 
clearly see a thin absorption feature, shaped like an incomplete 
elliptic ring encircling the waist of the nebula. It is observed in 
Ha (SuperCOSMOS image. Fig. 5 (a)\ in the near-IR (2MASS) 
and also at Spitzer 8.0 pm (Fig. 6) and 24 pm. The extinction is 
high on the lower side of the ring, probably corresponding to the 
front side (along the line of sight). This absorption is due to the 
dust present in the dense molecular material encircling the waist 
of the nebula. Two bright clumps, called Cl and C2 (Sect. 5.3; 
Fig. 7 (J)), are observed on each side at the waist of the neb¬ 
ula. We see, especially well at 5.8 pm or 8.0 pm (Fig. 6), several 
bright rims, probably bordering small size high density struc¬ 
tures embedded inside the two bright clumps observed at 250 pm 
(Fig. 5 (b)). A stellar cluster lies between these two condensa¬ 
tions, at the centre of the elliptic absorption ring. It has been 
detected on the 2MASS images by Dutra et al. (2003). Its centre 
lies at a(2000)= 15‘’03"'33*, ^(2000)= -57°40.1' (/=319?882, 
Z7=+ffi790). We discuss it in Sect. 5.2. 

The overall morphology of the G319.88+00.79 complex is 
discussed in Sect. 7.2. 


5.1. What more do we know about this region? 

It is a thermal radio source. The HI09a and HI 10a radio recom- 
binaison lines indicate a velocity of -38 km s“^ for the ionized 
gas (Caswell & Haynes 1987; the position observed, /=319?874 
does not correspond to the radio peak, but the large 
beam ~4'4 covers most of the radio emission); in the same di¬ 
rection V(H2C0)=-41 km s“^ (formaldehyde absorption). They 
conclude that G319.88+0.79 lies at the near kinematic distance 
of 2.6 kpc (using the Brand 1986 rotation curve). The near dis¬ 
tance is confirmed by the fact that this Hii region has an Ha 
counterpart; the distance is discussed in Sect. 5.2. The radio¬ 
continuum flux near 5 GHz is 1.6 Jy (Caswell & Haynes, 1987), 
and at 4.85 GHz it is 1.7 Jy (Kuchar & Clark 1997). Using 
Equation (1) in Simpson and Rubin (1990) and assuming an 
electron temperature of 5700 K as determined by Caswell & 
Haynes, we obtain an ionizing photon flux of 1.20x10^^ s“^; ac¬ 
cording to Martins et al. (2005), this points to an 08.5V exciting 
star, if it is single. The SUMSS map of the region shows a peak 
of radio emission in the direction of the central cluster, and a 
slight elongation in the direction of the two lobes (Fig. 5 (c)). 

The IRAS source IRAS 14597-5728 (coordinates 
/=319?8805, b=-\-(f.79ll) lies 5" away from the exciting 
star. Four IRDCs are listed by Peretto & Fuller (2009) in the di¬ 
rection or vicinity of G319.88+00.79 (Fig. C.l in Appendix C). 
They are discussed in Sect. 7.1. 
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Fig. 4. Bipolar nebula G319.88+00.79, com¬ 
posed of the S97 & S98 bubbles (Churchwell 
et al. 2006; the centres of the bubbles are in¬ 
dicated). (a) Composite colour image based on 
Spitzer observations. Red, green, and blue are 
for the 24 pm, 8.0 pm, and 4.5 pm emissions, 
respectively, (b) Composite colour image, with 
red, green, and blue for the 250 pm, 8.0 pm, 
and 4.5 pm emissions, respectively (logarith¬ 
mic units). The cold parental filament (or sheet) 
appears in red, and the PDR surrounding the 
ionized lobes in green 


5.2. The exciting duster - The distance of the complex 

A cluster is present in the central region (Fig. 5), visible in the 
near-IR on the 2MASS images and also in the Spitzer bands 
at 3.6 pm and 4.5 pm. The brightest star lies at /=319?8818, 
/?=-r0P7903; its 2MASS magnitudes are J = 10.476, H = 9.610, 
and K = 9.100. We have seen (Sect. 5.1) that the radio flux of 
the Hii region indicates an exciting star of spectral type 08.5V 
or more massive if ionizing photons are absorbed by dust. 
Assuming that this star is an 08.5V star, its near-IR photometry 
indicates that it is affected by a visual extinction of ~9.2 mag 
and lies at a distance of 2.1 kpc (using the photometry of O 
stars by Martins & Plez 2006 and the extinction law of Rieke & 
Lebofsky 1985). This confirms a near distance for this region; 
however, this distance is uncertain for two reasons: 1) ionizing 
photons can be absorbed by dust inside the ionized region (the 
24 pm emission comes from almost the same region than the 
radio continuum emission); if 70% of the ionizing photons are 
absorbed by dust we need an 07V exciting star. Thus a distance 


of 2.45 kpc. 2) Figure 5 shows {Inset) that the exciting star 
is double, composed of two stars separated by some 274; the 
2MASS magnitudes are thus overestimated. If the two stars were 
of equal intensity in K, their magnitude would be ^ = 9.85; this 
leads to a distance of 2.85 kpc. Thus in the following we keep 
the near kinematic distance of 2.6 kpc (+0.4 kpc)^. 


5.3. Molecular condensations 

Figure 1 (a) shows the dust temperature map of the region. The 
parental filament (or sheet) is cold with a minimum tempera¬ 
ture of 13.4 K in the direction of a region located along the fll- 


^The median kinematic distance uncertainty estimated for a large 
sample of H ii regions by Anderson et al. (2012) is 0.5 kpc. The distance 
uncertainty we estimate for this region is slightly smaller; the adopted 
distance range contains the two determinations of the photometric dis¬ 
tance. 
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Fig. 5. Centre of G319.88+00.79. (a) Colour image with - in red, green, 
and blue, respectively - the 4.5 fim emission, the K 2MASS emission, 
and the SuperCOSMOS Her emission of the ionized gas (linear units). 
(b) Colour image with 5.8 yum, 3.6 //m, and K 2MASS, in red, green, 
and blue, respectively (linear units). The 250 jum blue contours are pro¬ 
jected on the colour image (levels of 1000 to 7000 MJy/sr by steps of 
1000 MJy/sr). The exciting star is identified. The beam of the 250 jum 
image is at the lower left, (c) the radio continuum blue contours from 
the SUMSS survey at 843 MHz (emission of the ionized gas; levels 
of 0.010, 0.025, 0.050, 0.100, 0.200, and 0.300 Jy/beam; the SUMSS 
beam is at the lower left) are projected on the previous colour image. 
Inset 2MASS J band image of the exciting cluster (logarithmic units); 
the exciting star that appears elongated is probably double. 


ament. The dust in the PDRs surrounding the two ionized lobes 
is warmer with a maximum temperature of 23.2 K. 

When looking at the Herschel SPIRE maps, we see two 
bright clumps, Cl and C2. They are adjacent to the ionized re¬ 
gion, and their temperature is characteristic of PDRs (Anderson 




Fig. 6. Centre of G319.88+00.79 at 8.0 jjim showing the absorption fea¬ 
ture (the ellipse), the bright rims present at the waist of the nebula (green 
arrows), and the position of the exciting star. 


et al. 2011), respectively 22.3 K and 19.7 K (temperatures 
obtained by fitting the SEDs of these clumps, between 160 yum 
and 500 yum, with background subtraction). However, the 
column density map (Eig. 1(d)) shows the presence of another 
clump, C3, whose column density is as high as that of C2. 
Clump C3 is found at the coldest parts of the parental filament. 
It shows the dramatic influence of the dust temperature in 
the column density and mass estimates, hence the importance 
of the Herschel observations. The column density reaches 
values >8x10^^ cm“^ in the direction of C2 and C3. Clump 
C2, which is relatively warm (19.7 K), is clearly under the 
influence of the adjacent Hii region, whereas clump C3 (mean 
temperature 14.4 K) is a zone of the parental filament (or 
sheet) that lies farther away from the Hii region; it has pos¬ 
sibly not been influenced by it (as shown by its low temperature). 

The parameters of the three clumps are given in Table 2. 
The Galactic coordinates given in columns 2 and 3 are those 
of the peaks of column density. The peak’s column densities 
are given in column 4. The clumps Cl, C2, and C3 are the only 
structures of the field with a column density >5x10^^ cm“^. 
The dust temperatures, in Col. 5, are derived from the SEDs and 
based on the integrated fluxes of the clumps at 160 yum, 250 yum, 
350 yum, and 500 yum, corrected for the background emission. 
The uncertainty on T^ust given in Table 2 corresponds to an 
uncertainty of -15% on the fluxes^. To estimate the mass given 
in Col. 6 we have integrated the column density in an aperture 
corresponding to half the peak’s value. (These apertures are 
displayed in Eig. B.l, Appendix B.) We are thus dealing with the 
central regions of the clumps. The mass given in brackets has 
been estimated using the 250 yum flux integrated in an aperture 

^Note that the uncertainty on T^ust is possibly slightly higher. Using 
the dust opacity values of Ossenkopf & Henning (1994) leads to temper¬ 
atures that are higher by about 10%. The temperatures obtained without 
background correction are fainter by about 4% 
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Fig. 7. Temperature and column density maps of G319.88+00.79. The left side of the figure concerns the dust temperature, (a) The underlying 
grey image is the temperature map; the red contours correspond to temperatures of 19, 20, 21, and 22 K; the blue ones to temperatures of 14, 15, 
16, 17, and 18 K. The 36" beam of the temperature map (also that of the 500 jum map) is at the lower left, (b) The underlying image is the column 
density map in false colours, showing the cold parental filament, (c) The underlying image is the 8.0 fim emission map in false colours, showing 
the bipolar lobes. The right side of the figure shows the column density, (d) The underlying grey image is the column density map; the red contours 
correspond to column densities of 1, 2, 4, and 6 times 10^^ cm"^. The 18" beam of the column density map (also that of the 250 fim map) is at 
the lower left. The three clumps discussed in the text are identified, (e) The underlying image is the 250 jum emission map in false colours, (f) The 
underlying image is the temperature map in false colours. 


following the level at the 250 jim peak’s half-intensity and 
assuming the mean temperatures given in Col. 5. They differ by 
less than 20% of the masses estimated from the column density 
map. The uncertainty on the mass given in Col. 6 corresponds 
to the uncertainty on the temperature given in Col. 5^^. The 

^^The uncertainty in the mass is increased by the uncertainty on the 
distance. For example, a distance uncertainty of 0.4 kpc leads to a mass 


equivalent beam deconvolved radius at half-intensity is given in 
Col. 9. A mean density is estimated in Col. 10. (Using the mass 


uncertainty of ~75 M©. A larger source of uncertainty in the mass 
comes from the uncertain opacities. As an example, using the opaci¬ 
ties of Ossenkopf & Henning (1994) results in an increase in the dust 
temperature (24.6 K instead of 22.3 K for Cl) and therefore in a large 
decrease in the mass (150 M© instead of 246 M© for Cl). 
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and size of the central region and assuming an uniform density 
in this structure. This density is uncertain because we do not 
know the length of the clump along the line of sight and assume 
that the clump is spherical.) Clump C3 is far from spherical, and 
its size at half column density is 119.7" (along the filament) x 
24.3" (1.5 pc X 0.3 pc). If we assume that it is cylindrical (like 
a ID filament) with a length of 119.7" and a radius of 12.15", 
then its mean density is -8.0x10'^ cm“^. If its length along the 
line of sight is ~119.7" (more like a 2D structure), then its mean 
density is ~ 1.6x10^ cm“^. 


5 . 4 . Velocity field of the molecular clumps 

Two fields have been observed by MALT90 in this region: 
G319.901+00.792 and G319.872+00.788. They overlap and 
cover the Cl and C2 clumps, and the second field covers part 
of C3. (C3, less bright than Cl and C2 at 870 jum, has not been 
targeted by MALT90.) 

Figure 8 presents the radial velocity field of the Cl and C2 
clumps, obtained with MALT90 observations in the H^^CO^ (1- 
0) emission line. Similar results (not shown here) are obtained 
with the N 2 H^ molecular line. The velocity of the clumps lies in 
the range -41.5 to -44.5 km s“\ which is similar enough to that 
of the ionized gas (-38 km s“^), indicating that the molecular 
clumps and the Hii region are associated^Clump Cl and the 
filament (only a part of C3) have a mean velocity in the range 
-43.5 to -44 km s“^ and a line width in the range 3 to 4 km s“^ 
C2 has a mean velocity in the range -42 to -42.5 km s“^ 
Towards C2, the clump with the highest column density, the 
HCO^ line shows a double-peaked profile, which we believe is 
due to self-absorption. Nevertheless, because the velocity field 
shown is that obtained by doing a single Gaussian fit to a double- 
peaked profile, the resulting velocity lies close to the position of 
the self-absorption, so close to the expected systemic velocities 
of the clump (within 1 km/s). 

5 . 5 . Young stellar objects 

We used the colour-colour diagram [3.6]-[4.5] versus 
[5.8]-[8.0] displayed Fig. 9 to detect candidate Class I 
YSOs or some other noteworthy sources. The field considered 
here is centred on the central exciting star, and it has a radius 
of 7' to encompass the two lobes of the bipolar Hii region. 
The astrometry and photometry of these sources are given in 
Table 3, extracted from the Spitzer GLIMPSE Source Catalogue 
and completed by 24 pm and 70 pm measurements. The spectral 
index a and the colour [8]-[24] are also given in this table. 
Furthermore, a candidate Class 0 YSO is detected at 70 pm. All 
the sources discussed in the text are identified in Fig. 10. We 
first discuss two distinctive sources: 

• Source #3 appears as a stellar source on the Spitzer- 
GLIMPSE images, but as a region of extended emission on the 
MIPSGAL 24 pm image. Source 3 is one of the “intrinsically 
red” sources discussed by Robitaille et al. (2008), and according 
to these authors, it is very bright at 24 pm ([24] =0.129 mag). 
It is probably the emission of the extended region that has been 
measured. The 2MASS photometry of the central source points 
to an early B star (possibly B3V), affected by 5.8 mag of visual 


^^The velocity of the ionized gas has been obtained with a rather 
low resolution (~4'; Caswell & Haynes 1987), which prevents us from 
estimating the relative position of the clumps with respect to the ionized 
region. 



Fig. 8. Velocity field of the clumps located at the waist of 
G319.88+00.79. These results are obtained by MALT90 with the HCO^ 
(1-0) transition. What is shown here is the central velocity from doing a 
Gaussian fit to the line; only pixels with a ratio S/N>5 are shown here. 
The contours are for the column density (levels of 1, 1.5, 2, 3, 4.5, 7, 
and 8.5x10^^ cm"^). The 38" beam of the MALT90 observations is at 
the lower left. 
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[ 5 . 8 ] - [ 8 . 0 ] 


Fig.9. Colour-colour diagram, [3.6]-[4.5] versus [5.8]-[8.0], for all 
Spitzer-GhlMPSE sources with measurements in the four bands, and 
situated at less than 7' from the exciting star of G319.88+00.79. The 
big red dots correspond to sources with measurements more accurate 
than 0.2 mag in all four bands; the small red dots are for the less ac¬ 
curate sources. The blue dots are for candidate xAGB stars, follow¬ 
ing Robitaille et al. (2008). The location of YSOs dominated by a disk 
(Class II YSOs) and by an envelope (Class I YSOs) are indicated, fol¬ 
lowing Allen et al. (2004; see also Megeath et al. 2004). The extinction 
law is that of Indebetouw et al. (2005). The sources discussed in the 
text, mainly candidate Class I YSOs, are identified. 


extinction, if at the same distance than the bipolar nebula (an 
assumption). 

• Source #10 lies outside the top lobe of the bipolar nebula. 
It is most probably an evolved star, according to its 2MASS 
photometry, its brightness in the Spitzer-IRAC bands, and its 
low [5.8]-[8.0] and [8.0]-[24] colours. Also its spectral index 
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Fig. 10. Identification of the sources discussed in the text in the field 
of G319.88+00.79. (a) All sources located at less than 7' of the central 
exciting star (green plus). Source A is a compact 70 fim source. The 
underlying grey image is the Spitzer 5.8 pm image, (b) Sources located 
at the waist of the bipolar nebula and along the parental filament. The 
orange contours are for the column density (contours at 1, 2, 3, 4, and 5 
X 10^^ cm“^). The underlying grey image is the 70 pm image showing 
source A. 


is strongly negative (a < -2.5), and it has no detectable 70 pm 
counterpart. This evolved star is probably not associated with 
the bipolar nebula. 

Source #1 is located along the parental filament (on the left). 
It is probably a Class II YSO, according to its spectral index 
of -0.83 and colour [8]-[24]=2.36. It has no detectable 70 pm 
counterpart. Source #2 is a Class I YSO with a spectral in¬ 
dex of 1.58 and colour [8]-[24]=5.79; it is also detected and 
point-like at 70 pm and at 160 pm (flux of 0.77 Jy according to 
CuTEX). Source #2 is located in a direction of rather low col¬ 
umn density, 6.7x10^^ cm“^, thus a maximum external extinc¬ 
tion Ay ~7.2 mag leading to a corrected a ~1.4. This confirms 
that source #2 is a Class I YSO.) But the association of YSO #2 
with the bipolar nebula is highly uncertain because it lies outside 
the northern lobe. 

All the other candidate Class I YSOs, sources #4 to #9, lie 
along the parental filament. According to their colour [8] -[24] 
YSOs #4 and #9 are Class I, YSOs #5, #7, and #8 are flat- 
spectrum sources, and YSOs #6 is probably a highly reddened 
Class II. None of the YSOs #6 to #9 has a detectable 70 pm 
counterpart: it indicates that they are low-luminosity sources 
(<5Lo).YSO#4 lies near the centre of the nebula. It has a 70 pm 
counterpart, but its flux is very uncertain, indicating a luminosity 
~80 Lq. YSO #5 lies in the direction of the very centre of clump 
C2. Due to its position close to the PDR and close to bright rims 
emitting at 70 pm, we are unable to detect any counterpart. YSOs 
#6, #7, #8, and #9 lie on the border or centre of C3, in the direc¬ 
tion of the cold parts of the parental filament. It is also in this 
direction that a 24 pm and 70 pm compact source is detected 
without any counterpart at shorter wavelengths: we call it A. Its 
SED peaks at 70 pm and decreases at longer wavelengths. (Its 
160 pm flux is very uncertain, and it is not detectable at SPIRE 
wavelengths.) YSO A is probably the least evolved protostar. Its 


70 pm flux ~3.8 Jy indicates a central object with a luminosity 
~80 Lq. Eor this source, /6>g(/lE^(70)//lE^(24))=I.66±0.90, so it 
is possibly an extreme Class 0 YSO with a high infall rate (a 
PBR according to Stutz et al. 2013, Sect. 4.3). 

To conclude, except for YSO #2, which is possibly not as¬ 
sociated with the nebula, all candidate Class 0/I/flat-spectrum 
sources lie either at the waist of the nebula (YSOs #4 and #5) or 
along the parental filament. 


6. G010.32-00.15 

GOIO.32-00.15 is part of the well known bright star-forming re¬ 
gion W31, studied by Beuther et al. (2011) and presented as a su- 
perluminous star-forming complex. W31 contains numerous Hii 
regions. Two of them dominate the infrared and radio emission: 
GOlO.32-00.15 and GOlO.16-00.35 (see, for example, Kuchar 
& Clark 1997). GOlO.32-00.15 is clearly a bipolar Hii region. 
Eigure 11 (< 2 ) gives a Spitzer view of it; the 24 pm emission (red 
in the colour image) shows two lobes, the upper one open, the 
lower one possibly closed. The narrow waist is traced well by 
the bright 8.0 pm emission. This emission forms a closed ellipse 
(visible in Eig. 11 (b)), with a bright lower side and a faint up¬ 
per one. At this point it is difficult to know which one lies in the 
foreground and which in the background. 

GIO.32-00.15 is a radio-continuum source. The radio map 
at 21 cm obtained by Kim & Koo (2001) and displayed in 
Eig. \2 (a) clearly shows the bipolar structure of the ionized 
region. The radio emission resembles the 24 pm emission 
(Eig. 12 (b)) which is bright in the central region (where the 
24 pm emission is saturated), and which presents two diffuse 
elongated lobes. As discussed in Sect. 3, a fraction of the 24 pm 
emission comes from dust located inside the ionized region. 
The lobes seen at 70 pm (Eig. 12 (c)) surround the ionized 
lobes on the outside; this 70 pm emission is mainly due to dust 
located in the PDR surrounding the ionized region. The bipolar 
morphology of GO 10.32-00.15 has been discussed by Kim & 
Koo (2001, 2002). 

Molecular material is associated with this region. Kim & 
Koo (2002) mapped the region in (1-0) and CS (2-1) with 
a resolution of the order of 1'. Their map shows that: i) the 
molecular material associated with GOlO.32-00.15 and that as¬ 
sociated with GOIO. 16-00.35 share similar velocities, mainly in 
the range 10 to 13 km s“^ (their figure 5); ii) “the distribution of 
molecular gas is severely elongated in the direction orthogonal 
to the bipolar axis of GIO.3-0.1.” and their figure 9 illustrates 
this last point. Beuther et al. (2011) discuss and 

observations obtained with an angular resolution -27.5" and 
cold dust emission obtained at 870 pm (ATLASGAL survey) 
with a resolution -19". Their results agree with those of Kim & 
Koo for the velocity field. The high resolution of their 870 pm 
observations allowed them to detect four very bright cold dust 
clumps in the direction of GOlO.32-00.15 (their figure 1). These 
clumps are discussed in Sect. 6.3. 

Eigure \\ (b) presents a general view of the complex with 
the 250 pm emission of the cold dust in red, the 8.0 pm emission 
tracing the ionization front limiting the ionized material in 
green, and the stellar emission at 3.6 pm in blue. It shows the 
parental filament east of the nebula, and also the bright clumps 
present at the waist of the nebula. 
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Fig. 11. Bipolar nebula GOlO.32-00.15. (a) 
Composite colour image with red, green, and 
blue for the Spitzer 24 pm, 8.0 pm, and 4.5 pm 
emissions, respectively. (Note that the 24 pm 
emission is saturated in the central region), (b) 
Composite colour image with red, green, and 
blue for the 250 pm cold dust emission, 8.0 pm 
emission, and 3.6 pm stellar emission, respec¬ 
tively (the 250 pm emission is saturated in the 
C3 condensation, see Sect. 6.3). 


Seven IRDCs are listed by Peretto & Fuller (2009) in the 
vicinity of GOlO.32-00.15 (Fig. C.2 in Appendix C). They are 
discussed in Sect. 7.1. 

The overall morphology of the GOlO.32-00.15 complex is 
discussed in Sect. 7.2. 

6.1. Distance of G010.32-00.15 - Exciting cluster 

The distance of this region is very uncertain because different 
indicators give very different distances covering the range 


2 kpc to 19 kpc. This very peculiar situation is explained in 
Appendix A.l. In the following we try to identify the exciting 
star (or cluster) of GOlO.32-00.15 and determine its photomet¬ 
ric distance. 

What can we learn from the radio-continuum emission of 
the Hii region? The radio flux density of the Hit region is 
-14.7 Jy at 5 GHz (Du et al. 2011; Kuchar & Clark 1997). 
Assuming an electron temperature -7000 K (Reifenstein et al. 
1970), we estimate an ionizing photon flux (NLyc in s“^) such 
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Fig. 12. GOlO.32-00.15 at different wavelengths, showing the bipolar 
nature of the Hii region. The red contours correspond to the column 
density showing the dust clumps encircling the waist of the nebula (lev¬ 
els of 1.5, 1.0, 0.75, and 0.5x10^^ cm"^). The underlying grey images 
are (a) the VLA radio-continuum image at 21 cm, (b) the MIPSGAL 
map at 24 yum, and (cj the PACS map at 70 jum (all in logarithmic units to 
enhance the diffuse emission regions, especially the two ionized lobes). 


that log(NLyc)=48.65 or 49.72 for distances of 1.75 kpc or 6 kpc, 
respectively (justified in Appendix A.l), hence the spectral 
type of the exciting star, if single, respectively 07V or 03V 
according to Martins et al. (2005). For a far kinematic distance 
of 15 kpc the exciting cluster should contain more than seven 
03V stars. (Note that we do not take the possible absorption of 


ionizing photons by dust into account.) 
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Fig. 13. Exciting cluster of GOlO.32-00.15. We have identified the stars 
located at less than 10" of the exciting star (#1). (Their number in¬ 
creases with their distance to the exciting star.) The underlying grey 
image is the UKIDSS J image (logarithmic units). 


The deep near-IR UKIDSS images show that a cluster is 
present at the centre of the GOIO.32-00.15 Hii region (facing 
the C2 clump; Sect. 6.3). A star, our #1, dominates the cluster at 
all near-IR wavelengths (and up to 4.5 jum). Figure 13 displays 
the cluster in the J band. We have identified the 14 stars located 
at less than 10" of the exciting star, and brighter than 17.5 mag in 
J (for accurate photometry). Table A.l (Appendix A) gives the 
position of these stars, their JHK magnitudes according to the 
UKIDSS catalogue, and their distance d to star #1. Only the data 
for the central star (#1), saturated on the UKIDSS images, comes 
from the 2MASS catalogue or from Bik et al. (2005). Figure 14 
presents the (J - H) versus {H - K) and J versus (/ - K) dia¬ 
grams for all the stars located at less than 90" from the central 
star and for those at the centre of the cluster. The main sequence 
for a distance of 1.75 kpc, and a visual extinction of 18.4 mag is 
also drawn. (These distance and extinction are justified below.) 
Bik et al. (2005), from K-band spectroscopy, derive a spectral 
type 05V-06V for star #1. Using the synthetic photometry of O 
stars from Martins & Plez (2006) and assuming a spectral type 
05V for star #1, we derive a visual extinction of 18.8 mag and 
a distance of 1.6 kpc from the 2MASS magnitudes of this star 
(respectively 18.00 mag and 1.94 kpc from Bik et al. photome¬ 
try). This explains the adopted mean values used in Fig. 14. This 
figure shows that, except for two stars (#13 and #14), star #1 and 
surrounding stars are most probably main sequence stars, lying 
at the same distance thus forming a cluster; the visual extinction 
of the stars in the cluster is in the range 16.6 mag - 22.3 mag. 
Stars #13 and #14 are probably foreground stars. 

Star #1 dominates the cluster and is probably the main 
exciting star of the H ii region. Independent of its spectral type, 
its near-IR magnitudes can constrain the maximum possible 
distance of the cluster. Let us assume that star #1 is a 03V star; 
using 2MASS photometry, its visual extinction is 18.8 mag 
(independent of the distance), and its observed K magnitude 
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Fig. 14. Near IR photometry of the central exciting cluster of 
GOl0.32-00.15. (a) J - H versus H-K diagram of the cluster’s stars 
discussed in the text (photometry from the UKIDSS catalogue). The 
small red dots correspond to stars located at less than 90" of the ex¬ 
citing star (#1). The big red dots correspond to the stars identified in 
Fig. 13, located at less than 10" from the exciting star. The green plus 
correspond to O stars, affected by 0 and 18.4 mag of visual extinction. 
The blue line is the main sequence for B stars and stars of later spectral 
types. Main sequence stars are expected to lie in the region contained 
between the two black reddening lines (derived from main sequence O 
and MO stars), (b) J versus J-K diagram; the red symbols have the 
same meaning as in Fig. 14 (a). The two black lines limit the zone oc¬ 
cupied by reddened O stars at a distance of 1.75 kpc. The main sequence 
is also drawn for this distance, and for a visual extinction of 0 mag or 
18.4 mag (extinction law of Rieke & Lebofsky 1985). 


indicates a distance of 2.1 kpc (using the Bik et al.’s magnitudes, 
the maximum distance is 2.5 kpc). Since star #1 cannot be 
brighter than an 03V star, its distance cannot be greater than 
2.5 kpc. How can we reconcile a maximum distance of 2.5 kpc 
for GOlO.32-00.15 with the H 2 CO and Hi absorption spectra 


that point to greater distances? A tentative explanation is given 
in Appendix A. 1. 

To conclude, we have identified the exciting star of GO 10.32- 
00.14. It lies at the centre of a cluster, itself at the centre of 
the ionized region. This star is an 05V-06V star with a visual 
extinction of 18.4+0.4 mag. Its spectrophotometric distance is 
1.75+0.25 kpc. All this agrees with the radio-continuum flux 
(assuming that some ionizing photons are absorbed by dust) 
and with the near kinematic distance. We adopt this distance of 
1.75 kpc in the rest of this paper^^. 

6.2. Dust temperature and column density maps 

Figure 15 (< 2 ) presents the dust temperature map of this region. 
It shows that the warmer dust is found in the central region, 
close to the central ionized region, with the highest temperature 
(25.4 K) obtained in the direction of the PDR close to the central 
exciting cluster (at the edges of clumps C2 and C3, on their 
sides turned towards the exciting cluster and between them; 
the clumps are identified in Fig. 15 (d)). Low temperatures are 
found in the eastern filament (as low as 16.6 K), and in the 
direction of two clumps (as low as 16.2 K in clump C6 and 
16.6 Kin C5). 

Figure 15 (d) presents the column density map of this 
region. The column density is especially high over the region: 
higher than 2x10^^ cm“^ in the centre of the six clumps, Cl to 
C6 (Fig. 15 and Table 4). All these clumps (except C6) lie at the 
waist of the bipolar nebula. Clump C3 is saturated at 250 pm in 
its central region, and in unsaturated pixels the column density 
reaches 3.8x10^^ cm“^. A hole in the column density is present 
in the central region (in the direction of the exciting cluster) 
where the column density is of the order of 2x10^^ cm“^. The 
ionized gas lies in this central region, so this column density 
is probably that of the foreground or background extended 
emission. (We adopt this value for subtracting the background 
to estimate the parameters of the clumps in Sect. 6.3.) 

In this region we again see the influence of the temperature 
on the determination of the column density. Clumps C5 and C6 
are not among the brightest at 250 pm, but they are cold, and as 
a consequence their column density is high. 


6.3. Molecular clumps 

The field of Fig. 15 contains eight structures with a peak col¬ 
umn density higher than 1.5x10^^ cm“^. The five brightest are 
located at the waist of the nebula; we call them Cl to C5. Cl 
and C3 are almost circular, and Cl is the most compact. All 
are adjacent to the ionization front and bordered by a bright 
rim at 8.0 pm (on the side turned towards the central exciting 
cluster). The parameters of the high column density clumps can 
be found in Table 4. The first line of Cols. 2 and 3 gives the 
Galactic coordinates of their column density peaks. The second 
line of Cols. 2 and 3 gives the coordinates of the ATLASGAL 
emission peaks at 870 pm^^, while Col. 4 gives their velocity, as 

^^At this near distance, W31 is no longer one of the most lumi¬ 
nous giant molecular cloud/high-mass star formation complexes in the 
Galaxy, as claimed by Beuther et al. (2011). 

^^The positions of the column density peaks and of the ATLASGAL 
870 pm emission peaks given by Beuther et al. differ by less than 4". 
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Fig. 15. Dust temperature and column density maps of GOlO.32-00.15. The three left images deal with the temperature; the green cross shows the 
position of the exciting star, (a) The underlying grey image is the temperature map; the red contours correspond to temperature of 21, 22, 23, and 
24 K, the blue ones to temperatures of 19, 18, and 17 K, the yellow one to 20 K. (b) The same contours are superimposed on the Spitzer 8.0 pm 
image, which enhances the waist of the nebula, (c) The same contours are superimposed on the column density map. (d) Column density map 
(contour levels of 0.3, 0.5, 0.75, 1.0, 1.5, and 2.0 xlO^^ cm"^). The clumps discussed in the text are identified. 


obtained from the NH 3 emission (Wienen et al. 2012) or from 
the emission (Beuther et al. 2011). Column 5 gives their 
peak column density. Column 6 gives - only for information - 
the mean dust temperature estimated from the temperature map 
(Fig. 15). The mass, given in Col. 7, has been obtained by inte- 


This is compatible with the uncertainties on positions expected for these 
two sets of observations obtained with beams ~19". 


gration of the column density in an aperture following the level 
corresponding to half of the peak value (after correction for the 
background). This mass is the mass of the central region of the 
clump. Column 8 gives the size of the central region, in arc- 
sec (and parsec in brackets), estimated as the equivalent radius 
of the apertures used (hence the beam-deconvolved size at half 
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Fig. 16. Identification of sources discussed in the text, in the vicinity 
of GOlO.32-00.15. (a) Column density map showing the five clumps 
present at the waist of the bipolar nebula (red contour levels at 0.3, 0.5, 
0.75, 1.0, 1.5, and 2.0 xlO^^ cm"^). (b) Spitzer 5.8 pm image; the plus 
symbols show the positions of the methanol masers, of the EGOs, the 
MYSOs, the UC Hii region, and IRAS 18060-2005; the exciting star 
is also identified. The orange arrows show the pressure exerted by the 
ionized gas trying to escape from the central ionized region through 
the molecular shell, between the dense clumps, (c) MAGPIS image at 
21 cm showing the UC Hii region and the dense ionized layers bor¬ 
dering the IF (especially the borders of C2 and C3 turned towards the 
exciting cluster). 


peak value). Column 9 gives an estimate of the mean density of 
the central region, assuming that it is homogeneous. 

The central regions of the C2 to C4 clumps are rather similar 
in terms of size, with a radius in the range 0.10 pc - 0.12 pc, 
a peak column density a few 10 ^^ cm“^, a mass in the range 
220 Mq - 330 Mo, and a mean density of several 10^ cm“^. 
They are dense. The central region of Cl is not very different, 
but it is smaller and denser, so it deserves to be called a core. 
Clumps C5 and C 6 are also dense and are slightly more extended 


and more massive than the Cl to C4 clumps. Clump Cl which 
is not associated with the region (on the basis of its velocity, 
see below), lies at an unknown distance, so we cannot estimate 
its mass and density. Two velocity components are observed in 
the direction of the Cl extension (Sect. 6.4); thus, two structures 
are present along this line of sight, and we do not know which 
fraction of the molecular material is associated with the bipolar 
nebula. 

The Cl to C5 clumps form a structure that limits, at least on 
the southern side, the brightest part of the ionized region. Their 
velocities are similar to the velocity of the ionized gas, as shown 
by Table 4 and discussed in Sect. 6.4. Figure 16 (Jb) shows that 
the IF is distorted between the clumps. We suggest that this 
is due to the high-pressure ionized gas trying to escape from 
the central region via low-density paths between the clumps. 
All clumps display signatures of massive-star formation at 
work (see Sect. 6.5); each of the Cl to C4 clumps contains a 
6.67 GHz methanol maser in its very centre, at less than 3" 
(0.025 pc) from the column density peak (Green et al. 2010). 
Two of these masers in Cl and C4 lie close to an extended green 
object (EGO; Cyganowski et al. 2009), another one in C3 lies 
close to an UC Hii region (Walsh et al. 1998), and the last one 
in C2 lies close to a massive YSO (YSO #4). Another methanol 
maser lies in the direction of the Cl extension. All these objects 
are identified in Fig. 16 and are discussed in Sect. 6.5. 


6.4. Velocity field of the molecular clumps 

Seven overlapping fields have been observed by MALT90 that 
cover this bright molecular complex. Three main velocity com¬ 
ponents are present along the line of sight in different zones of 
the field. 

• The first component, corresponding to velocities in the 
range 10.5 km s“^ to 13.5 km s“^ as shown in Fig. 17, is ob¬ 
served in the direction of the clumps Cl to C5 (those at the waist 
of the bipolar nebula) and in the direction of C 6 . These veloc¬ 
ities are similar to that of the ionized gas, showing that these 
clumps and the central Hii region are parts of the same com¬ 
plex. It was already known, but indirectly, that C 6 is in contact 
with the bipolar nebula, because the southern lobe is distorted 
at the border of C 6 . These velocities are also in good agreement 
with the velocities measured with the NH 3 line by Wienen et 
al. (2012) and with CO lines by Beuther et al. (2011, and refer¬ 
ences therein). Self-absorption is clearly observed in the direc¬ 
tion of C3 (Fig. 17 Right), the clump with the saturated 250 pm 
emission; these lines present a double-peaked profile over the C3 
clump, the result of a central self-absorption. The H^^CO^ line, 
which is optically thin, shows only one component centred on 
the self-absorption feature in the direction of C3 (Fig. 17 Right). 

• The second component, corresponding to velocities in the 
range 28 km s“^ to 32.5 km s“^ as shown in Fig. 18, is observed 
in the direction of Cl and vicinity. 

• The third component, corresponding to velocities in the 
range 43.5 km s“^ to 47 km s“^ as shown in Fig. 19, is observed 
east of Cl and C2 (especially in the direction of the Cl exten¬ 
sion structure), and south of C4, in directions adjacent to these 
clumps. 

The results presented in Figs. 17, 18, and 19 have been 
obtained with the H^^CO^ (1-0) or HNC (1-0) molecular lines; 
similar results are obtained with these two lines. (The velocities 
agree to within 0.5 km s“^) 
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Fig. 17. The main velocity component observed over the GOlO.32-00.15 field. Left, velocity field obtained with the HCO^ (1-0) molecular line. 
The contours are for the column density (levels of 0.5, 0.75, 1, 1.5, and 2x10^^ cm"^). The 38" beam is at the lower left. Right velocity profile 
from a region of 36"x36"centred on C3, at /=10P299 The HCO^ line shows self-absorption whereas the line is optically thin. 

The red dotted line shows the velocity obtained by fitting the splitted HCO^ line with one gaussian. 


These observations show that at least three molecular clouds 
are present in the field at velocities of the order of 12 km s“\ 
30 km s“\ and 45 km s“^ thus presumably at different distances. 
These clouds are observed here via their emission, and are prob¬ 
ably the same clouds that are responsible for the H 2 CO absorp¬ 
tion of the radio-continuum emission of the central H ii region 
(Appendix A.l). 
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Fig. 18. Second velocity component observed over the GOlO.32-00.15 
field (results obtained with the HNC (1-0) molecular line). The contours 
are for the column density (levels of 0.4, 0.5, 0.75, and 1x10^^ cm"^). 
The black crosses indicate the position of the MAGPIS UC H ii region. 
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Fig. 19. Third velocity component observed over the GOlO.32-00.15 
field (results obtained with the HCO^ (1-0) molecular line). 


6.5. Young stellar objects 

In the following we try to detect candidate Class O/I YSOs 
and to confirm their evolutionary status. Figure 20 presents the 
Spitzer [3.6]-[4.5] versus [5.8]-[8.0] diagram: the red dots are 
for sources located less than 7' from the exciting star #1, thus 
in a field covering the bipolar nebula, lobes included, its associ¬ 
ated molecular clumps, and the parental filament. The blue dots 
correspond to candidate xAGB stars. GOlO.32-00.15 lies in the 
direction of the galactic bulge, so it is not surprising to find many 
candidate AGB stars in this field. We have identified in Fig. 20 
all the sources that we consider as candidate Class I YSOs or that 
we discuss in the following. Their parameters, name (column 1), 
coordinates (columns 2 & 3), magnitudes from the near-IR to 
70 pm (columns 4 to 12), luminosity (column 13), spectral in¬ 
dex a between K and 24 pm (column 14), and colour [8]-[24]^^ 

note a discrepancy between the 24 pm magnitudes in the 
Robitaille et al. (2008) and the MIPSGAL catalogues (see Table 5). The 
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Fig.20. Colour-colour diagram, [3.6]-[4.5] versus [5.8]-[8.0], for all 
Spitzer-GIAMPSE sources with measurements in the four bands, situ¬ 
ated at less than 7' from the exciting star (#1) of GOlO.32-00.15. The 
big red dots correspond to sources with measurements more accurate 
than 0.15 mag in all four bands, the small red dots for the less accurate 
sources. See Fig. 9 for further information. 






(column 15) are given in Table 5. YSOs #16 to #26 are located 
far from the waist of the bipolar nebula, and many of them are 
possibly not associated with it. Several candidate Class I YSOs 
are located in the vicinity of the nebula, especially in the di¬ 
rection of the clumps and filaments present at the waist of the 
nebula. We discuss the content of each of the clumps in the fol¬ 
lowing. (In this section the column densities are given without 
background subtraction.) 

• Cl: Figure 21 shows clump Cl and its vicinity. YSO #1, 
with [8]-[24] in the range 4.16 mag to 4.49 mag, is a Class I 
YSO. It is located outside of Cl but on the upper side of the waist 
of the bipolar nebula, in a region of intermediate column density 
(3.3x10^^ cm"^), and has a faint 70 pm counterpart (indicating 
a luminosity ~12 L©). 

Figure 2\ (b) Sc (d) show that several sources lie at the very 
centre of Cl (at less than 3" - 0.025 pc - from the column den¬ 
sity peak): 1) an EGO (Cyganowski et al. 2008, 2009), which 
appears as a rather bright green source in Fig. 21 (a)\\i presents 
a 24 pm counterpart; 2) a class II methanol maser (Walsh et al. 
1998, Green et al. 2010), almost coincident with the EGO. Its 
velocity, 15.4 km s“^ (Green et al. 2010), indicates that it is as¬ 
sociated with Cl. A class I methanol maser is also associated 
with the EGO and the class II methanol maser (Cyganowski et 
al 2009). A bright and compact 70 pm source lies in the very 
centre of Cl. It is the third-brightest source of the field at this 
wavelength with a fiux density of 197+20 Jy. Its luminosity, es¬ 
timated from its fiux at 70 pm, is -1500 L©. It lies in the exact 
direction of the EGO and the methanol maser and has no radio 
counterpart. We suggest that Cl contains a MYSO in an early 
evolutionary stage because it has no near-IR counterpart (thus 


consequence is that some sources classified (according to their [8]-[24] 
colour) as flat-spectrum sources using the Robitaille et al. magnitudes 
are classified as Class II sources when using the MIPSGAL magnitudes. 


a massive Class O/I YSO). Two slightly extended 8.0 pm and 
24 pm sources (sources A & C in Fig. 21 (c) are also present in 
the direction of Cl. The bright source A lies -375 - 0.3 pc - from 
the centre of Cl. But the 70 pm emission peaks in the direction 
of the EGO & methanol masers rather than in the direction of 
source A. Sources A & C have no detectable radio counterpart 
(fiux <5 mJy on the MAGPIS map at 20 cm). 

Cl has an extension of high column density. In this direction 
there are two YSOs: YSO #2, with a colour [8]-[24]=3.55 mag 
and a faint 70 pm counterpart is a fiat-spectrum YSO. YSO #3, 
with a colour [8]-[24]=5.03 mag, and a bright 70 pm counter¬ 
part is clearly a Class I YSO. A 6.7 GHz methanol maser has 
been detected in its direction by Green et al (2010), which lies at 
less than 1" of YSO #3. YSO #3 is slightly extended at 4.5 pm, 
suggesting that it can be associated with an EGO. However YSO 
#3 is most probably not at the same distance as the bipolar 
GOlO.32-0.15 Hit region. Its associated methanol maser has a 
velocity of 49.9 km s“^ showing that YSO #3 is associated with 
the molecular material of the third velocity component (near 
44 km s“^) that overlaps with the Cl extension (see Sect. 6.3.1). 
A slightly extended source, source B in Fig. 21 (h), (c), (d) is 
also present in the direction of the Cl extension. (It lies at 372 
from the column density peak.) It is seen at 8.0 pm, 24 pm, and 
70 pm and has a very faint radio counterpart on the MAGPIS 
map at 20 cm, so it is probably a faint UC H ii region. We esti¬ 
mate its fiux at 20 cm to be ~8 mJy, and it is unresolved. Its fiux 
density at 70 pm is -38 Jy. We do not know if it is associated 
with the bipolar nebula or with the molecular component at a 
velocity -44 km s“^ (same situation for YSOs #2). 

• C2: C2 faces the exciting cluster. Fig. 23 (a). It is bor¬ 
dered by a bright emission region on the side facing the cluster. 
It is bright at radio wavelengths (MAGPIS image; Fig. 23 (c)), 
where we see the emission of the dense ionized layer border¬ 
ing the molecular clump, and it is bright at 8.0 pm or 5.8 pm 
(Fig. 23 (b)), since we see the emission of the PAHs in the PDR, 
excited by the radiation leaking from the H ii region. 

A 6.67 GHz methanol maser lies at the very centre of C2 
(Walsh et al. 1998; Green et al. 2010). It lies at 274 - 0.2 pc - 
from the column density peak. Its velocity, 11.5 km s“^ (Green 
et al. 2010), indicates that it is associated with the clump. In 
the same direction is YSO #4, a very bright near- and mid-IR 
unresolved source. YSO #4 is not part of the GLIMPSE cata¬ 
log, probably because it lies in a region of bright nebular emis¬ 
sion and is saturated in several GLIMPSE bands. (This region is 
also saturated at 24 pm.) YSO #4 is the second brightest 70 pm 
source of the field (after the UC H ii region in C3). It is point-like, 
with a fiux of 630+50 Jy, indicating a luminosity -4600 L©. It 
has no detectable radio counterpart on the MAGPIS image (fiux 
<5 mJy at 20-cm). Source #4 lies in a massive condensation 
(-400 M©) and shows methanol maser emission at 6.7 GHz. All 
this indicates that source #4 is a massive YSO (MYSO), which 
possibly is accreting material too strongly to form a detectable 
Hii region (to explain the lack of an associated UC Hii region; 
Walmsley 1995). It has a near-IR counterpart, so it is probably a 
massive Class I and not a Class 0 YSO. 

The nature of YSO #5 is not clear. It lies in the direction of 
the central exciting cluster, 12" from star #1 (Fig. 22). Its JHK 
magnitudes correspond to a B3V star with about 16.8 mag of 
visual extinction, thus possibly part of the exciting cluster. Its 
8.0 pm counterpart is bright and slightly extended. YSO #5 lies 
in the saturated 24 pm zone. As shown by Whitney et al. (2013), 
an early B star illuminating an optically thin dust cloud can ap¬ 
pear to have an IR excess and resemble an YSO. It is possibly 
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Fig. 21. Cl clump and surroundings. The 
sources discussed in the text are identified, (a) 
Colour image with the 5.8 jjm, 4.5 jim, and 
3.6 fim emissions in red, green, and blue, re¬ 
spectively (logarithmic units). The red contours 
refer to the column density (levels of 0.3, 0.5, 
0.75, 1.0, 1.5, and 2 10^^ cm"^). (b) Colour im¬ 
age with the 70 //m, 8.0 jum, and 3.6 jum emis¬ 
sions in red, green, and blue, respectively (lin¬ 
ear units), (c) MIPSGAL image at 24 //m; the 
contours are for the 70 fim emission (levels at 
5000, 10000, 15000, and 20000 MJy/sr). A, B, 
and C are slightly extended sources of 8.0 yum, 
24 yum, and 70 jum emission, (d) Herschel 
70 jim image. 


the case for source #5, superimposed on a relatively bright back¬ 
ground at 8.0 /xm. 



Fig. 22. Identification of YSO #5. The underlying grey image is the 
Spitzer 8.0 pm image. The contours are for the column density (levels 
of 0.5, 1.0, and 1.5x10^^ cm“^); they show the centre of clump C2. 
YSO #4 is also identified (saturated with two glints). 

• C3: C3 is bordered by a bright rim, which is especially 
bright on its east side facing the exciting cluster (Fig. 23 (a)). 
Bright radio-continuum emission is observed in the direction of 
this bright rim. Figure 23 (c) also shows an UC H ii region in the 
direction of the centre of C3. We suggest that it is a separate, in¬ 
dependent, HII region embedded inside C3 or lying on the bor¬ 


der of C3. Its size is 0.06 pc (FWHM, beam deconvolved). Its 
flux is given in Table A.2 (Appendix A). Its ionizing photon flux 
(log(NLyc)~47.2) corresponds, according to Smith et al. (2002), 
to a main exciting star (if single) of spectral type BOV - B0.5 V. A 
small cluster is present in this direction, visible on the UKIDSS 
K image and on the Spitzer images up to 5.8 pm, which is prob¬ 
ably the exciting cluster. But we have not been able to identify 
the exciting star of this UC H ii region. This UC H ii region has 
a mid-IR counterpart, a small extended region of similar size 
emitting at all GLIMPSE wavelengths. (This region is saturated 
at 24 pm.) It is the brightest 70 pm compact source of the fleld 
with a flux of 890+30 Jy. A faint class II methanol maser (Walsh 
et al. 1998; Green et al. 2010) has been detected in the direction 
of C3. It lies at 1.874 - 0.016 pc - from the column density peak 
and at 675 from the centre of the UC Hii region. It is probably 
associated with C3 since its velocity is 19.9 km s“^. 

Another small region of extended GLIMPSE emission lies 
in the direction of C3 (region D, Fig. 23). This one has no de¬ 
tectable radio-continuum counterpart (flux <5 mJy), but it has a 
saturated 24 pm counterpart and a faint extended 70 pm counter¬ 
part. Again a small cluster is present in its direction. It contains 
a star that could be a late B star with a visual extinction of about 
22 mag (UKIDSS catalog: /= 17.089 mag, 14.795 mag, 
K=13.354 mag; coordinates /=10?2941 b=-(f.l443). This re¬ 
gion is probably located in front of C3 because, in its direction, 
the column density is of the order of 2x10^^ cm“^, which corre¬ 
sponds to a visual extinction of about 215 mag, so much higher 
than that of the star. 

Two faint candidate Class I YSOs (#6 and #7) lie outside of 
C3 but still in a region of high column density and along a 8.0 pm 
filament originating in the waist of the bipolar nebula. The col- 
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umn densities in the direction of #6 and #7 are 5.2x10^^ cm“^ 
and 9.7x10^^ cm“^, respectively. YSO #6 has no measurable 
counterparts at 24 yum (due to saturation) and 70 yum; its nature 
is uncertain. YSO #7 has a faint 24 yum counterpart and no de¬ 
tectable one at 70 yum. Its nature is uncertain (a low-luminosity 
Class Il/flat-spectrum YSO). 




Fig. 23. Clumps C2 and C3. (a) Colour image with the Spitzer 5.8, 4.5, 
and 3.6 pm emissions in red, green, and blue, respectively (logarithmic 
units). The red contours are the column density contours (same levels as 
in Fig. 21(b)). The candidate Class I YSOs #4, #5, #6, and #7 are iden¬ 
tified. The positions of the methanol masers are indicated by plusses. 
(b) Colour image with the 70 pm, 5.8 pm, and 3.6 pm respectively in 
red, green and blue (linear units), (c) The red column density contours 
are superimposed on the MAGPIS radio-continuum image. The UC H ii 
region is identified. 


• C4 and C5: Fig. 24 shows clumps C4, C5, and vicinity. 
In the very centre of C4 (at the peak of column density) lies 
an EGO (Cyganowski et al. 2008), which is much fainter than 


the one in Cl, and a class II methanol maser brighter than the 
one in Cl (Walsh et al. 1998; Green et al. 2010). The velocity 
of the maser, 4.6 km s“\ indicates that it is associated with C4. 
It lies at r .'8 - 0.015 pc - from the column density peak. The 
EGO and the class II methanol maser have no Spitzer counterpart 
(except at 4.5 pm\ but no 24 pm emission), but they have a faint 
70 pm counterpart, as shown by Eig. 24 (b). C4 with its EGO 
and methanol maser, and with no near-IR and radio counterparts, 
resembles Cl. However, its 70 pm emission is not as bright as the 
Cl one. (Its flux of 30+10 Jy indicates a luminosity ~260 L©.) 
It is a less luminous Class O/I YSO than the one in Cl. 

C5 lies west of C4. Two sources lie in the direction of C5: 
1) YSO #8, a flat-spectrum or a Class II YSO with [8]-[24] in 
the range 2.95 mag - 3.29 mag, which has no detectable 70 pm 
counterpart; 2) a small source of extended emission, which we 
call E, visible at all GLIMPSE wavelengths and at 24 yum and 
70 pm (Eig. 24). Source E resembles source B, with a 70 pm flux 
~24 Jy, but it has no detectable 20-cm counterpart (flux <5 mJy). 
We suggest that it is a small region of PAH and thermal emission 
around a late B star, but this star is not detected in UKIDSS or 
Spitzer-GLlMPSE images. The column density in the direction 
of YSO #8 and of region E is in the range 1.8 - 2.7x10^^ cm“^. 

A faint 8.0 pm filament extends to the south of C4 and C5, 
originating in the waist of the bipolar nebula. Two YSOs, #9 and 
#10, are located along this filament. The column density in the 
direction of this filament and of YSOs #9 and #10 is in the range 
9 - 10x10^^ cm“^. YSO #9, with [8]-[24] in the range 2.62 mag 
- 3.36 mag is a Class II or a flat-spectrum YSO. It has a faint 
70 pm counterpart, indicating a luminosity ~50Lo. YSO #10 
is faint at all Spitzer wavelengths and has no detectable 70 pm 
counterpart, so it is possibly a low-luminosity Class II YSO. Two 
velocity components, near 13 km s“^ and 46 km s“\ are present 
in the direction of YSO #9 and #10, making their association 
with the bipolar nebula uncertain. 

Two more clumps lie close to the bipolar nebula at the border 
of the southern lobe: 

• C6: it is a cold clump that belongs to the complex. Its veloc¬ 
ity is similar to the velocities of the ionized gas and of clumps Cl 
to C5. C6 is similar to the clumps Cl to C5 in terms of column 
density, mass, and density. C6 contains three candidate Class I 
YSOs, #11 in its very centre (at 270 - 0.017 pc - from the column 
density peak), and #12 and #13 nearby (Eig. 25). YSO #11 is 
bright at 24 pm (with [8.0]-[24]=6.06 mag) and bright at 70 pm, 
so it is clearly a Class I YSO. Its 70 pm flux ~30 Jy indicates a 
luminosity -260 L©. YSO #12 is almost as bright as YSO #11 at 
GLIMPSE wavelengths. YSO #13 is brighter than #11 and #12 
from K to 4.5 pm, but much fainter at 8.0 pm. Both YSO #12 
and #13 lie in the direction of the diffraction ring of YSO #11, 
so their 24 pm magnitude is very uncertain. YSO #12 is possi¬ 
bly a flat-spectrum source. The nature of YSO #13 is unknown. 
No radio emission is detected on the MAGPIS radio map in the 
direction of C6 (flux <5 mJy at 20-cm). 

C6 was probably a pre-existing cloud, reached by the IP dur¬ 
ing the expansion of the southern lobe (as C6 protudes inside the 
southern lobe). It interacts with the southern lobe and distorts its 
edge. C6 appears to be compressed on its north-east border by 
the ionized gas. Pigure 25 shows that the molecular material is 
not symmetrically distributed around YSO #11: we measure a 
column density gradient -6.3x10^^ cm“^ per parsec on the side 
of C6 turned towards the ionized gas and a column density gra¬ 
dient of -2.6x10^^ cm“^ per parsec on its side turned towards 
the outside. 

• C7: this clump does not belong to the GOlO.32-00.15 
complex because its velocity is too different from that of the 


21 










L. Deharveng et al.: Bipolar nebulae 




- 00 . 10 “ 


(U 

■O -00.11“ 


ro -00.12“ 

03 


-00.13“ 


Galactic Longitude 


-00.14“ 

10.31“ 


10.30“ 


10.29“ 


10.28“ 


10.27“ 


Fig. 24. C4 and C5 clumps and vicinity, (a) Same colours and units as 
in Fig. 21 (a). Sources #8, #9, and #10 are candidate Class I YSOs. The 
brightness has been enhanced in the vicinity of the EGO. The slightly 
extended GLIMPSE source E is identified, (b) Same colours, units, and 
contours as in Eig. 21(b). The positions of the EGO and of the methanol 
maser are indicated. 


ionized gas of the bipolar Hii region. The candidate Class I 
YSOs found in the direction of C7 and the UC Hii regions 
present there are discussed in Appendix A.2. 

Several other candidate Class I YSOs lie in the vicinity of 
GOIO.32-00.15. They are identified in Appendix A, Fig. A.2. 
Four of them lie in the direction of the parental filament east 
of the bipolar nebula. They are most probably part of the com¬ 
plex. The column density is high in their direction (Fig. 26) with 
N(H2)>7 x10 ^^ cm“^. YSO #20 has a bright 24 jim counter¬ 
part ([8]-[24] in the range 2.87 mag to 3.36 mag) and a faint 
70 yum one. It is a low-luminosity Class II or flat-spectrum YSO 
(L~6Lo). At 70 yum, YSO #20 has a “twin” of similar bright¬ 
ness, YSO #21, that is very faint at 24 yum. These two sources 
are separated by about 16". YSO #21 is possibly a Class 0 YSO 
of low-luminosity (~9 L©). But it is not an extreme Class 0 
(a PBR according to Stutz et al. (2013). YSO #22 has a faint 
24 yum counterpart and no detectable 70 yum one. Its colour, 
[8]-[24]=2.36 mag, indicates that it is a Class II YSO. YSO #23 
is a bright 24 yum and 70 yum source. Its colour, [8]-[24] in 
the range 4.50 mag to 4.83 mag, and its 70 yum flux density of 



10.25“ 10.24“ 10.23“ 10.22“ 10.21“ 10.20“ 

Galactic Longitude 


Fig. 25. C6 dust clump and vicinity. Colour composite image with the 
24 yum, 8.0 yum, and 4.5 yum emissions in red, green, and blue, re¬ 
spectively (linear units). The three candidate Class I YSOs are iden¬ 
tified. The red contours are for the column density (same levels as in 
Fig. 21 (b)). 


10.5 Jy indicate a Class I of intermediate luminosity (~100 L©). 
Furthermore, YSO #23 has bright 160 yum and 250 yum point-like 
counterparts (with fluxes of 20.1 Jy and 18.5 Jy, respectively); 
we estimate the mass of its envelope to be 17 M© (for a temper¬ 
ature of 18.2 K taken from the temperature map). The YSOs #20, 
#21, #22, and #23 are most probably associated with the bipo¬ 
lar nebula. But they lie far from the IF limiting the H ii region 
(the closest, YSO #23, is located in projection some 1.4 pc away 
from the IF), so their formation has probably not been triggered 
by the central nebula. 

The other candidate Class I YSOs are not observed in the 
direction of structures linked to the nebula, and their velocity 
is unknown. Their association with GOlO.32-00.15 is therefore 
uncertain. They are discussed in Appendix A.2. 

To summarize, among the sources selected as candidate 
Class I YSOs on the basis of their Spitzer-IRAC colours and 
clearly associated with the bipolar H ii region, we And four con¬ 
firmed Class I YSOs (#1, #4, #11, and #23), three flat-spectrum 
or Class II sources (thus of uncertain nature between Class I 
and Class II; #8, #12, and #20), and one Class II YSO (#22). 
Furthermore, YSO #21, in the direction of the parental filament, 
is detected as a Class 0, and the clumps Cl and C4 contain 
Class 0 sources detected thanks to their associated outflows (via 
EGOs and class II methanol masers). All the high-luminosity 
sources (bright at 70 pm) are observed in the direction of the 
massive clumps forming the waist of the nebula. 

Several small-size diffuse regions, regions A to F, are present 
in the vicinity of the clumps, generally in the direction of 
rather high column density regions. They are detected at Spitzer- 
GLIMPSE wavelengths, especially at 8.0 pm. They have coun¬ 
terparts at 24 pm and 70 pm. We suggest that they are associated 
with B stars, mostly stars later than B2 because they have no 
detectable radio-continuum emission (region B is an exception), 
but stars that are hot enough to heat the very small dust grains 
that emit at 24 pm and possibly at 70 pm. 
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Fig. 26. YSOs in the direction of the eastern parental filament. The red 
contours correspond to column densities of 3, 4, 5, 6, and 7 10^^ cm"^. 
They are superimposed on grey images (logarithmic units): (a) the 
24 jum image; (Z?) the 70 jum Hi-GAL image. The 18" beam of the col¬ 
umn density map (also that of the 250 jum image) is at the lower left. 


7. Discussion 

7.1. IRDCs and star formation 

IRDCs were first identified as regions of low brightness at mid- 
IR wavelengths, this property being attributed to absorption of 
the background emission by dense dusty clouds (Perault et al. 
1996; Egan et al. 1998; Carey et al. 1998). Infrared dark clouds 
(IRDCs) are often considered as the birthplace of massive stars 
or clusters (for example, Battersby et al. 2010). Thousands of 
IRDCs have been detected in the Galactic plane (Simon et al. 
2006; Peretto & Fuller 2009). 

When considering 3171 IRDCs that were ''Spitzer-dsirk re¬ 
gions”, Wilcock et al. (2012) have shown that only 1205 of them 
were ''Herschel-hvight clouds”, hence dense structures seen in 
emission at SPIRE wavelengths and not simply “holes” in the 
sky, zones devoid of material. This mis-identification of IRDCs 
occurs in our regions; for example, four IRDCs are catalogued 
in the G319.88-h00.79 field, and three of them do not corre¬ 
spond to dense structures, as shown by Fig. C.l (Appendix 
C). They correspond to the centres of the ionized bubbles de¬ 
void of 8.0 pm emission (SDC G319.882-h0.810 at the cen¬ 
tre of S98 and SDC G319.886-h0.753 at the centre of S97), 
or to a zone of low 8.0 pm emission between two filaments 
(SDC G319.859-h0.818). Only one of them corresponds to a 
zone of real absorption owing to the presence of a dense structure 
along the line of sight (SDC G319.910-h0.796; on the left border 
of Cl). Seven IRDCs are listed in the vicinity of GOlO.32-00.15, 
as shown by Fig. C.2 (Appendix C). Five of them are observed 
in the direction of high column density structures, especially 
SDC GOlO.215-0.206 in the direction of the C6 clump and 
SDC GOlO.358-0.166 -h SDC GOlO.405-0.200 in the direction 
of the parental filament. But two of them are not at all associ¬ 
ated with dense structures: SDC G010.315-0.133 which lies in 
the direction of the central hole in column density (a zone oc¬ 


cupied by the ionized region and devoid of 8.0 pm emission); 
SDC GOlO.291-0.187 which lies in the direction of the centre 
of the bottom lobe (also a zone occupied by ionized gas). 

The very dense clumps that lie at the waist of 
GOIO.32-00.15 and contain massive sources presently forming 
or just formed are not detected as IRDCs. This situation is 
frequently observed (see also the C2 clump in G319.88-f00.79, 
or the massive-star forming condensation 1 in RCW120, 
Deharveng et al. 2009, Zavagno et al. 2010). These massive 
clumps, not detected as IRDCs and containing MYSOs, are 
adjacent to the PDRs of Hit regions, and the absorption they 
induce is hidden at near or mid-IR wavelengths by the bright 
emission of the nearby PDRs. 

7.2. The global morphology of G319.88-\-00.79 and 
G010.32-00.15 - Filaments and sheets 

Many studies are devoted to filaments and to star formation 
along filaments. It is often forgotten that 2D structures - sheets, 
viewed edge-on, can be mistaken as filaments. Also perfect 
geometrical structures (cylinders or flat planes) are never met 
in nature where density gradients and wavy structures are 
everywhere present. What have we learned about the overall 
morphology of the G319.88-f00.79 and GOlO.32-00.15 bipolar 
nebulae and their parental clouds? 

That we are seeing a ring (or torus) of dense absorbing ma¬ 
terial at the waist of the G319.88-h00.79 bipolar H ii region, cen¬ 
tred on the exciting star(s) (Sect.5; Fig. 6), indicates two things: 
1) the exciting star(s) formed in a 2D dense structure, a sheet 
rather than a cylindrical filament; 2) this parental sheet is not 
seen edge-on; the line of sight makes an angle of about 58° with 
the symmetry axis of the bipolar nebula (assuming a circular 
ring; the observed size of the elliptic ring is ~ 1.95' x 1.03', thus 
1.47 pc X 0.78 pc for a distance of 2.6 kpc). Figure 27 shows 
our view of the morphology of this region. Only the bottom part 
of the ring is clearly seen in absorption, which implies that this 
absorption is due to material located in front of the H ii region. 
We do not detect Herschel-SPIRE emission associated with this 
structure, which suggests that there is little material. The bright 
rims observed left of the ring are associated with the Cl clump, 
and the bright rims observed right of the ring are associated with 
C2. We suggest that C2 lies slightly at the back of the H ii region 
and that it is composed of several small-size dense cores (each 
bordered by the bright rims observed, for example, at 8.0 pm. 
Fig. 6), not resolved by Herschel-SPIRE. Thus dense structures 
are located at the back of the H ii region, but we do not know the 
extent of the molecular material there. 

The mean velocity of the ionized region is -38 km s“^ 
(Caswell & Haynes 1987), and the velocity of the molecu¬ 
lar material is in the range -44 km s“^ to -42 km s“^ (Cl 
to C3 clumps). If the morphology illustrated by Fig. 27 is 
correct, then the ionized bottom lobe is dominant at radio 
wavelengths. (In the lobes the ionized gas flows away from the 
dense central region.) Observations of the velocity field of the 
ionized region is necessary to confirm the proposed morphology. 

In the case of the GOIO.32-00.15 bipolar Hii region we 
know a little more about the velocity field of the ionized gas. 
The velocity of the ionized gas has been measured by Araya et 
al. (2007) at three positions with a beam of 2.5'. In the direc¬ 
tion of the central region (/=10?313, b=-(f.l52), the velocity is 
9.4 km s“^; the velocity is 12.6 km s“^ 2.5' north in the direction 
of the upper lobe; the velocity is 5.6 km s“^ 2.5' south of the 
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Fig. 27. Morphology of G319.88+00.79. The ionized material appears 
in pink, the neutral one in blue. We do not know the extent of the molec¬ 
ular material at the back of the nebula. 


central region in the direction of the bottom lobe. This velocity 
gradient indicates an approching bottom lobe with respect to the 
central region, and a receding upper one (thus Fig. 28). 

The parental filament (left of the H ii region) is seen in ab¬ 
sorption, and in line with it, the Cl and the C4 + C5 external 
clumps are also seen in absorption at Spitzer-\KAC wavelengths. 
Thus these clumps are slightly in front of the H ii region. The C2 
and C3 clumps are not clearly seen in absorption, but they are 
bordered by bright rims (PAH emission) on their side turned to¬ 
wards the exciting star (Fig. 16); thus they lie in front of the 
ionized region and its exciting star^^. This suggests that the Cl 
to C5 clumps form a single structure, half a torus of dense mate¬ 
rial, encircling the waist of the bipolar H ii region and located in 
front of this region. 

The H II region is not open on its back side (it is ionization- 
bounded), as we see a faint Spitzer-IKAC PAHs emission closing 
the clumpy torus on its back side. The parental structure is much 
more dense on the front side of the H ii region than on the back 
one. We suggest that the parental structure was a fiat 2D cloud 
presenting a strong density gradient or containing a dense pre¬ 
existing filament. 

7.3. Triggered star formation 

Several processes can trigger star formation in the vicinity of 
HII regions, as illustrated by Deharveng et al. (2010; their figure 
4): small-scale gravitational instabilities of the IF, large-scale 
gravitational instabilities of the material accumulated in a shell 
behind the IF (the so-called “collect and collapse process”), and 
radiation-driven implosion (“RDI process”) of a pre-existing 
dense clump compressed by the ionized gas. The process 
proposed by Fukuda & Hanawa (2000; Sect. 3) is somewhat 
different, involving the gravitational collapse of clumps, formed 
by the compression of a pre-existing dense filament due to 
an adjacent Hii region. It is often difficult to establish which 
triggering process is at work. And, as shown by the simulations 
of Dale et al. (2012; 2015), it is often difficult to deduce 
which objects have been induced to form and which formed 
spontaneously. 


^^If they were lying behind the exciting star, their whole area would 
show a bright emission at 8 pm. Here we see only their top border turned 
towards the exciting star. 



Fig. 28. Morphology of the GOlO.32-00.15 complex. The plane repre¬ 
sented here contains the parental filament, the clumps at the waist of the 
nebula, and the exciting star. The line of sight makes a small angle with 
this plane. 


We propose the GOIO.32-00.15 region as an example of 
massive-star formation triggered by an Hii region and involv¬ 
ing different processes. Concerning the Cl to C5 clumps, the 
signatures of triggering follow. 

• An incomplete shell of dense molecular material surrounds 
the ionized region at its waist. It is fragmented, and the five 
fragments are massive (several hundreds of solar masses). The 
shape of the PDR shows the interaction of the high-pressure ion¬ 
ized gas with the surrounding molecular material: the ionized 
gas tries to escape from the central H ii region via low density 
zones present between the dense fragments (Fig. 16). 

• Each of the fragments contains at least one embedded 
young stellar object, most often of high mass. Cl contains a 
MYSO in an early stage of evolution (probably a Class 0), C2 
contains a high-luminosity YSO (probably a Class I) with a 
high accretion rate, C3 contains an UC H ii region, and C4 con¬ 
tains a YSO probably of Class 0. All these massive YSOs have 
methanol maser emission at 6.7 GHz. According to Breen et al. 
(2010), this maser emission indicates an age <5x10^ yr for the 
YSOs. We cannot determine the dynamical age of the bipolar 
H II region accurately because appropriate models are missing, 
but we can estimate a lower limit for its age. The waist of the H ii 
region has a diameter of 1.44 pc (maximum size in the parental 
structure). An Hii region of radius 0.72 pc, which is excited by 
an 05V - 06V star (NLyc=10^^), and is evolving in a homo¬ 
geneous medium of density 10"^ H 2 cm“^ (a reasonable density 
for the parental structure) has an age -170000 yr (according to 
the evolution of an H ii region described by Dyson & Williams 
1997). A bipolar Hii region should expand more slowly in the 
plane of the parental structure than a spherical Hii region. A 
fraction of the ionized gas fiows away from the dense central 
region towards the lobes, reducing the pressure of the ionized 
gas in the central region. The bipolar H ii region is most proba¬ 
bly older than 170000 yr, hence much older than the Class O/I 
YSOs located in the clumps at its waist. 

• The Cl to C5 clumps form a single structure: a dense half 
shell encircling the waist of the nebula. The probability that these 
clumps were pre-existing structures distributed along a ring and 
simultaneously reached by the expanding IF is very low. We ar¬ 
gue that their formation results from the interaction of the H ii 
region with the dense surrounding material. Under the assump¬ 
tion that triggering is indeed taking place, there are two alter- 
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native solutions for the triggering process. The first option is 
the equivalent of the collect and collapse process, but taking 
place in a 2D rather than in a 3D cloud. The material of the 
parental sheet is collected during the expansion of the central 
Hii region, forming a dense torus. Later, gravitational instabili¬ 
ties lead to the fragmentation of the torus, followed by star for¬ 
mation. However, at present, no theoretical model or simulation 
can replicate this scenario. The second option is that, according 
to Fukuda & Hanawa’s similations, the H ii region compressed a 
nearby pre-existing filament, leading to the formation of the C2 
and C3 central cores that contain the more evolved and massive 
YSOs (YSO #4 and an UC Hii region). The external clumps. Cl 
and C4 -h C5, contain the less evolved YSOs, which are detected 
via the presence of EGOs, methanol masers, and 70 jim emis¬ 
sion. The present data do not allow us to determine that Cl and 
C4 -r C5 are first-generation cores, or if only C2 and C3 are first- 
generation cores while Cl and C4 -h C5 are second-generation 
ones. However, some observational details do not agree with the 
simulations. In Fukuda & Hanawa’s simulations, the filament is 
not distorted by the H ii region. (It does not encircle the waist 
of the Hii region.) In the simulations, the Hii region is den¬ 
sity bounded on the side opposite to the filament, whereas our 
HII regions are ionization bounded in the plane of their parental 
structure. 

The C6 clump suggests that another triggering process 
may be at work. It is clearly a massive pre-existing structure 
that has been reached by the IF and probably compressed by 
the high-pressure ionized gas. An intermediate-mass Class I 
YSO possibly results from this compression, which is located 
at the column density peak of the clump. We estimate that 
its luminosity is -260 L©. We suggest that its formation re¬ 
sults from the radiation driven implosion of the pre-existing 
dense cloud, as described and simulated by numerous authors 
(Bisbas et al. 2009, 2011, and references therein). Bisbas et al. 
(2011) have shown that the incident ionizing fiux is the critical 
parameter determining a cloud evolution. The surface of C6, 
located in projection at about 2.9 pc from the 05 - 06 exciting 
star, receives a ionizing photon fiux -10^ cm“^ s“^ which is 
a case for low photon fiux. The shock front progressing in the 
cloud is weak, and it sweeps up matter slowly. According to 
the simulations, star formation occurs in a dense filament, well 
behind the IF, towards the centre of the cloud. The observations 
of C6 fit this picture well; however, the resolution of the 
Herschel observations does not allow access to the small-scale 
structure of the clump. Furthermore, the clouds involved in the 
simulations of Bisbas et al. have much lower mass than the C6 
clump (-600 Mq, obtained by integration in a circular aperture 
of radius 50"covering the whole clump). Thus triggered star 
formation in C6 via the RDI process also needs confirmation. 

Star formation is less active around G319.88-h00.79 than 
around GOlO.32-00.15 (Sect. 5.4). We suggest that the forma¬ 
tion of YSOs #4 and #5, located at the waist of the nebula, near 
the ionization front, could - based on their location - have been 
triggered by the central H ii region. It is especially the case of 
YSO #5, which is located in the direction of a bright rim (part of 
the C2 clump), thus close to the IF. All the other YSOs in this re¬ 
gion are located along the parental filament, outside the Cl and 
C2 clumps. The Class 0 YSO A is located about 0.6 pc (in pro¬ 
jection) from the IF at the right waist of the nebula. The same 
situation is observed in the GOlO.32-00.15 field; for example, 
YSOs #21 and #23 are located along the parental filament, too 
far away from the IF for triggering. (The Class I YSO #23 is 
located about 1.4 pc away - in projection - from the IF at the 


left waist of the nebula; the Class 0 YSO #21 lies even farther 
away.) The formation of these Class O/I YSOs has most likely 
not been triggered by the expansion of the bipolar nebulae. Their 
formation along a filament probably results from some other pro¬ 
cess, such as the filament’s fragmentation via gravitational in¬ 
stabilities (Jackson et al. 2010; Andre et al. 2010 and references 
therein). 


8. Conclusions 

Using the Spitzer GLIMPSE and MIPSGAL surveys and the 
Herschel Hi-GAL survey we have identified 16 candidate bipo¬ 
lar nebulae along the Galactic plane, between -h60° and -60° 
in Galactic longitude, -rl° and -1° in latitude. This first paper 
focuses on two of them, G319.88-1-00.79 and GOIO.32-00.15, 
which are especially interesting for illustrating the morphology 
of these regions and the formation of stars in their vicinity. Our 
conclusions follow. 

• Concerning their morphology both regions present similar 
signatures. 

- The parental structure, whether filament or sheet, is cold 
(13 K-17 K), and dense (-10"^ cm“^ or more in G319.88-r00.79). 

- The brightest parts of the H ii region and its exciting cluster 
are centred on the parental structure. 

- Two ionized lobes, closed or open, extend perpendicular to 
the parental structure. Warm dust (20 K-25 K) is present in the 
PDRs bordering the ionized regions, including the lobes. 

- High column density clumps are present at the waist of 
the nebulae, adjacent to the ionized regions. MALT90 obser¬ 
vations show that they are associated with the ionized regions, 
because they present similar velocities. Their mass is a few hun¬ 
dred solar masses. High densities are measured in their central 
regions (several 10"^ cm“^ in G319.88-h00.79, several 10^ cm“^ 
in GOlO.32-00.15). 

• We identified the exciting stars of the two Hii regions. 
This has been especially useful for estimating the distance of 
GOlO.32-00.15. The distance of this region, as adopted by dif¬ 
ferent authors, is in the range 2 kpc to 19 kpc. We have estimated 
a spectrophotometric distance of 1.75+0.25 kpc, which is com¬ 
patible with the near kinematic distance. 

• We used colour-colour diagrams created from the Spitzer- 
GLIMPSE observations to identify candidate Class I YSOs in 
the vicinity of the Hii regions. Some of these have been con¬ 
firmed, using their 24 pm and Hi-GAL 70 pm emissions. Class 0 
YSO were detected thanks to their 70 pm emission. 

Near G319.88+00.79, eight candidate Class I YSOs were 
identified by their Spitzer-GIAMPSE colours. Thanks to their 
[8]-[24] colour, three of them (YSO #2, #4, and #9) are con¬ 
firmed Class I, three are fiat-spectrum sources (thus of uncertain 
evolutionary status between Classes I and II; YSOs #5, #7, #8), 
and two are probable Class II (YSOs #1 and #6). A Class 0 is 
detected (YSO A), which has a luminosity -SOL© and proba¬ 
bly a high infall rate. Only two YSOs (#4 and #5) are located 
in a clump at the waist of the nebula. The other YSOs, except 
YSO #2, are located along the parental filament. 

Near GOIO.32-00.15 twenty-six sources are candidate 
Class I YSOs, according to their Spitzer-GLIMPSE colours; 
among them, eight are confirmed Class I YSOs (#1, #3, #4, #11, 
#15, #23, #24, and #25), five are Class II YSOs, and ten have 
uncertain natures. Class II or fiat-spectrum sources, due to un¬ 
certain 24 pm magnitudes. Three Class 0 are present (2 of them, 
in Cl and C4, with associated EGOs and methanol masers and 
without near-IR or Spitzer-IRAC counterparts, except at 4.5 pm; 
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and YSO #21). Several of these YSOs are possibly not associ¬ 
ated with the bipolar nebula. This region demonstrates the im¬ 
portance of velocity measurements to associate a clump or a 
YSO with a nearby Hii region. All the high-luminosity sources, 
bright at 70 yum, are observed in the direction of the massive 
clumps forming the waist of the nebula. The two brightest and 
possibly most evolved sources (the UC H ii region in C3, and the 
Class I YSO #4 in C2, of luminosity -4600 L©) are embedded in 
the central C2 and C3 clumps. The less luminous sources, also 
possibly the least evolved. Class 0 YSOs (associated with EGOs 
and class II methanol masers), are located in the outside Cl and 
C4 clumps. 

• The GOlO.32-00.15 bipolar Hii region has triggered the 
formation of stars in the surrounding medium via various pro¬ 
cesses. We suggest that the C6 massive clump was a pre-existing 
cloud reached by the ionization front limiting the southern lobe 
and compressed by the high-pressure ionized gas. An interme¬ 
diate luminosity Class I YSO (-260 L©) is located at its centre, 
probably formed by triggering via the RDI process. A cluster 
is possibly forming there because two other YSOs are located 
nearby. 

On the contrary, we suggest that the Cl to C5 massive 
clumps were not pre-existing clouds, since they form a single 
circular structure, but that their formation results from the in¬ 
teraction of the Hii region with the surrounding dense molec¬ 
ular material. This triggering can result from the collect of the 
dense material of a flat cloud (thus the formation of a dense 
torus at the waist of the H ii region) or via the compression of 
a pre-existing dense filament as suggested by Fukuda & Hanawa 
(2000). However, no model or simulation is presently available 
to confirm one or the other process. 

• GOlO.32-00.15 is more extreme than G319.88-h00.79. Its 
exciting star is more massive (05V-06V versus 08V-09V), 
and the condensations present at its waist are more numerous 
(5 against 2), more compact (equivalent radius in the range 
0.05 pc - 0.18 pc against -0.30 pc), and denser (in the range 
3x10^ cm“^ - 5x10^ cm“^ against less than 10^ cm“^). More 
Class O/I YSOs are present in the vicinity of GOlO.32-00.14 
than near G319.88-h00.79, and several of them are massive. The 
high density of the molecular material into which the exciting 
star of GOlO.32-00.15 formed and evolved^^ could be the 
reason for such a situation, but only a statistical approach can 
confirm or disprove this assumption. 

After this work was completed, we became aware of the pa¬ 
per by Dewangan et al. (2015), a study of star formation around 
the GOlO.3-0.1 Hii region. These authors base their estimate 
of the evolutionary status of the sources on their near-IR and 
Spitzer-GLIMPSE magnitudes. They do not mention that there 
are several velocity components in some directions, hence the 
uncertain association of many sources with the central Hii re¬ 
gion. Furthermore, and this is an important point, they mention 
the bipolar nature of the central H ii region, but in their discus¬ 
sion of (triggered) star formation they do not take the very spe¬ 
cific morphology of this star-forming complex into account. 
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Table 2. Properties of the clumps in the field of G319.88+00.79. 


Name 

1 

b 

N(H2)' 

T ^ 

^ dust 

m2 

Req" 

n(H2)2 


n 

n 

(cm-2) 

(K) 

(Mo) 

(") 

(cm"^) 

C 1 

319.9012 

+00.7935 

5.3x1022 

22.3+1.9 

303 (246;y2) 

25.2 (0.32 pc) 

5.4x10“* 

C2 

319.8668 

+00.7868 

8.3x1022 

19.7+1.4 

497 (458;f2o) 

24.9 (0.31 pc) 

8.6x10“* 

C3 

319.8324 

+00.7925 

7.5x1022 

14.4+1.2 

766 (623;™) 

119.7x24.3 (1.5 pcxO.3 pc) 

1.6-8.0x10“* 


^ Peak column density corrected for a background value ~4.6xl0^^ cm“^ (the lowest value of column density in a field of radius 10' centred on the exciting 
star; this low column density value is located at /=320°04 b=+(f.lS, outside any feature associated with the bipolar nebula and its parental cloud). 

^ Dust temperature derived from the SED, based on the integrated fluxes of the clumps in the range 160 yum-500 /im. 

^ Mass obtained by integrating the column density in an aperture following the level at the N(H 2 ) peak’s half intensity. In brackets is the mass obtained by 
integrating the 250 /rm emission in an aperture following the level at the 250 /rm peak’s half intensity and using the temperature in column 5. 

Equivalent radius of the aperture following the level at the N(H 2 ) peak’s half intensity (beam deconvolved). 

^ Mean density of the central region, enclosed by the aperture following the level at the N(H 2 ) peak’s half intensity. 


Table 3. Sources in the field of G319.88+00.79, mainly candidate Class I YSOs or objects discussed in the text. 


Name 

1 

o 

b 

o 

J 

(mag) 

H 

(mag) 

K 

(mag) 

[3.6] 

(mag) 

[4.5] 

(mag) 

[5.8] 

(mag) 

[8.0] 

(mag) 

[24] 

(mag) 

S(70yum) 

(Jy) 

a 

[8]-[24] 

(mag) 

1 

319.95786 

+00.80312 

14.334+0.060 

12.813+0.035 

11.980+99.999 

10.879+0.053 

10.109+0.056 

9.551+0.037 

8.992+0.033 

6.63 ROB 

<0.2* 

-0.83 

2.36 

2 

319.91510 

+00.87341 

15.521+0.478 

14.278+0.358 

13.042+0.137 

11.532+0.086 

10.679+0.044 

9.063+ 0.030 

7.225+0.020 

1.43 ROB 

2.97+0.19 CuTEx 

1.58 

5.79 

3 

319.89360 

+00.76700 

12.029+0.024 

11.430+0.027 

11.134+0.021 

10.901+0.047 

10.706+0.068 

10.274+0.102 

8.967+0.129 

(0.129 ROB) 

<0.3* 



4 

319.87887 

+00.78600 




10.961+0.156 

9.701+0.199 

8.378+0.065 

7.584+0.079 

3.76+0.3* 

4+2* 


3.82 

5 

319.86893 

+00.78639 

17.314+99.999 

16.758+99.999 

13.791+0.048 

10.457+0.167 

9.249+0.071 

7.694+0.092 

6.96+0.10* 

3.54+0.3* 


1.06 

3.42 

6 

319.85309 

+00.78660 




12.578+0.058 

11.545+0.079 

10.732+0.081 

10.416+0.139 

>8.0* 

<0.2* 


<2.4 

7 

319.84955 

+00.78840 




12.859+0.096 

10.988+0.069 

10.006+0.055 

9.301+0.047 

5.71 ROB 

<0.2* 


3.59 

8 

319.84122 

+00.79177 




14.011+0.092 

12.724+0.115 

11.891+0.131 

11.502+0.109 

8.00+0.50* 

<0.2* 


3.50 

9 

319.83133 

+00.78281 




13.286+0.064 

11.528+0.060 

10.391+0.058 

9.590+0.028 

5.13+0.15* 

<0.2* 


4.46 

10 

319.78050 

+00.84672 

8.057+0.024 

6.459+0.017 

5.749+0.016 

7.207+0.278 

5.641+0.066 

5.286+0.033 

5.231+0.021 

4.92+0.15* 

<0.2* 

-2.57 

0.31 

A 

319.85210* 

+00.78850* 








5.98+0.15* 

3.88+0.24 CuTEx 




Notes. Our own measurements are indicated by an *. ROB indicates 24 /urn magnitudes from Robitaille et al. (2008). CuTEx indicates values 
taken from the Hi-GAL catalogue, to be published (Molinari et al. 2015, in prep.). They have been obtained using the package CuTEx (Molinari 
et al. 2011); for an example of using CuTEx, see Elia et al. (2013). 
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Table 4. Molecular clumps in the field of GOlO.32-00.15 (C7 and a fraction of the Cl extension are not associated 
with the bipolar H ii region). 


Name 

/> 


Vlsr 

N(H2)‘ 

T/ 

Mass^ 

Req" 

n(H2) 




(kms-‘) 

(10^^ cm"^) 

(K) 

(Me) 

(") (pc) 

(10^ cm"^) 

Cl 

10.3430 

-00.1418 


3.27 

19.6 

184 

5.7 (0.05) 

55 


10.3419 

-00.1422 

12.14? 10.1’ 






C2 

10.3232 

-00.1600 


2.09 

22.0 

223 

14.4 (0.12) 

4.2 


10.3224 

-00.1604 

12.02? 11.9’ 






03^ 

10.2998 

-00.1460 


>2.66 

21.8 

330 

12.6 (0.11) 

9.2 


10.3001 

-00.1469 

12.83? 11.2’ 






C4 

10.2868 

-00.1249 


2.80 

19.4 

216 

11.7 (0.10) 

7.6 


10.2876 

-00.1243 

11.7’ 






C5 

10.2837 

-00.1137 


2.50 

17.1 

506 

21.0 (0.18) 

3.1 


10.2836 

-00.11.36 

10.6’ 






C6 

10.2282 

-00.2067 


1.76 

16.3 

250 

18.7 (0.16) 

2.2 


10.23 

-00.21 

12.14® 






Cl 

10.3215 

-00.2572 


1.61 

19.4 


9.9 



10.32 

-00.26 

32.15® 






Cl extension 

10.3551 

-00.1470 


1.57 

19.0 


19.9 



^ The coordinates (first line)and the column densities are those of the column density peaks measured on the column density map. The peak 
column densities have been corrected for a background emission (~0.2xl0^^ cm“^, measured in the central region near /=l(k318 Z?=-0?146). 
The coordinates (second line) are those of the ATLASGAL 870 /rm emission peaks from Beuther et al. (2011). 

^ Mean temperature measured using the dust temperature map by integration in the aperture corresponding to the central region. 

^ Mass of the central region of the clumps obtained by integrating the column density in an aperture following the level at the N(H 2 ) peak’s half 
intensity. 

Equivalent radius of the aperture following the level at the N(H 2 ) peak’s half intensity (beam deconvolved). 

^ This clump is saturated at 250 /rm; the coordinates and column densities have been measured on the column density map based on the 350 /rm 
emission map; due to a smoothing by the beam (larger at 350 /rm than at 250 /rm), the peak column density is underestimated. 

^ V(NH 3 ) from Wienen et al. (2012); pointed observations in the direction of the ATLASGAL 870 jjm intensity peaks; FWHM 40". 

^ V(C^^O) from Beuther et al. (2011); spectra extracted from a velocity map towards the positions of the ATLASGAL 870 /rm intensity 
peaks; FWHM 27!5. 
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Table 5. Candidate Class I YSOs and other sources discussed in the field of GOlO.32-00.15. 


Name 

1 

b 

J 

H 

K 

[3.6] 

[4.5] 

[5.8] 

[8.0] 

[24] 

S(70jum) 

L 

a 

[8]-[24] 


(°) 

(°) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(Jy) 

{Lo) 


(mag) 

Cl 

1 

10.34286 

-0.12324 


16.906+0.044 

14.980±0.016 

12.330±0.070 

10.764±0.050 

9.066±0.052 

7.199±0.024 

3.03 (2.71 ROB) 

1.1±0.3* 

11.6 

1.84 

4.17 (4.49 ROB) 

2 

10.35481 

-0.15608 



15.932±0.038 

11.921±0.116 

11.096±0.125 

8.136±0.062 

6.469±0.114 

2.92 

6.2±2.0* 


2.21 

3.55 

3 

10.35612 

-0.14847 




11.930±0.087 

9.595±0.079 

8.721 ±0.044 

8.962±0.162 

3.93±0.10* 

41±10* 



5.03 

EGO 

10.3423 

-0.1422 








1.81±0.20* 





methanol maser 

10.3422 

-0.1422 













methanol maser 

10.3561 

-0.1483 













A 

10.3439 

-0.1415 













B 

10.3547 

-0.1463 









38±10* 




C 

C2 

4 

10.3460 

-0.1432 













10.3226 

-0.1596 


16.55* 

13.08* 

7.86* 

saturated 

saturated 

saturated 

saturated 

630±50* 

4600 



5 

10.31474 

-0.15225 

14.220±0.002 

12.447±0.001 

11.424±0.001 

10.429±0.077 

9.403 ±0.145 

7.599±0.113 

5.846±0.217 

saturated 





methanol maser 
C3 

6 

10.3226 

-0.1597 













10.29520 

-0.15858 




12.125+0.147 

10.293±0.081 

8.945±0.075 

8.461±0.174 


<1.0* 

<10.6 



7 

10.28971 

-0.16641 



17.500+0.159 

13.440+0.096 

11.022±0.05 

9.169±0.051 

7.894±0.032 

4.5±0.5* 

<0.5* 

<5.5 


3.4 

methanol maser 

10.2993 

-0.1461 













UC HII region 

10.3008 

-0.1471 









890±30* 




D 

10.2941 

-0.1443 













C4 and C5 















8 

10.28264 

-0.11271 



14.772+0.013 

9.864+0.039 

8.573±0.038 

7.357±0.03 

6.609±0.026 

3.66 (3.32 ROB) 

<1.0* 

<10.6 

1.52 

2.95 (3.29 ROB) 

9 

10.27116 

-0.12776 




12.221±0.061 

9.599±0.054 

7.759±0.039 

6.920±0.033 

4.30 (3.56 ROB) 

5.1±1.0* 

49 


2.62 (3.36 ROB) 

10 

10.26919 

-0.13257 




14.014±0.137 

11.980±0.074 

10.252±0.084 

8.910±0.055 

6.0±0.5* 

<4.5* 



2.9 

EGO 

10.2874 

-0.1247 













methanol maser 

10.2872 

-0.1252 













E 

C6 

11 

10.2873 

-0.1088 









24±5* 




10.22781 

-0.20707 



15.225±0.020 

12.261 ±0.067 

10.841 ±0.079 

8.744±0.034 

7.021 ±0.022 

0.96 

30.0±5.0* 

262 

2.64 

6.06 

12 

10.22664 

-0.20910 



16.572±0.067 

12.473 ±0.069 

10.234±0.052 

8.799±0.041 

7.684±0.035 

4.2±1.0* 



1.65 

3.5 

13 
Cl 

14 

10.23044 

-0.20471 


17.006+0.049 

13.151+0.003 

10.678±0.049 

9.443 ±0.042 

8.865±0.039 

8.569±0.063 






10.32254 

-0.22372 




13.702±0.122 

12.032±0.141 

10.558±0.088 

10.070±0.070 

>6.7* 

<0.4* 



<3.4 

UC HII region 

10.3202 

-00.2331 









130±10* 




15 

10.32035 

-0.25867 




11.412±0.119 

8.849±0.113 

7.231±0.090 

6.002±0.296 

1.71 

113±30* 



4.30 

methanol maser 

10.3204 

-0.2585 













HII region 

10.3218 

-0.2587 













E 

10.3151 

-0.2547 









14±3* 




External 















16 

10.25432 

-0.24416 

18.969±0.121 

16.088±0.021 

13.083±0.003 

10.163±0.070 

9.114±0.073 

8.234±0.038 

7.706±0.036 

4.56±0.10* 

1.6±0.2* 


0.39 

3.15 

17 

10.27684 

-0.23287 

19.278±0.161 


15.234±0.020 

12.800±0.069 

11.641 ±0.074 

10.697±0.092 

9.865±0.057 

>7.5* 

<0.3* 



<2.4 

18 

10.36044 

-0.23713 



15.595±0.027 

11.595 ±0.060 

10.249±0.047 

9.217±0.045 

8.137±0.032 

5.61 (5.31 ROB) 

0.2±0.1* 


1.07 

2.53 (2.83 ROB) 

19 

10.35452 

-0.22285 

18.489±0.079 

16.901±0.045 

15.479±0.025 

12.443 ±0.062 

11.182±0.076 

10.188±0.055 

9.453±0.051 

6.44 ROB 

<0.3* 


0.59 

3.01 ROB 

20 

10.40812 

-0.20145 



16.939±0.095 

12.015±0.086 

9.839±0.058 

8.549±0.039 

7.633±0.027 

4.76 (4.27 ROB) 

0.6±0.2* 

6.5 

1.99 

2.87 (3.36 ROB) 

21 

10.4038 

-0.2028 








7.47 ROB 

0.8±0.3* 

8.6 



22 

10.38881 

-0.19460 


16.443±0.029 

13.718+0.005 

11.327±0.051 

10.315±0.050 

9.590±0.048 

9.261 ±0.074 

6.9±0.5* 

<0.3* 

<3 

-0.26 

2.36 

23 

10.37284 

-0.18482 



14.631+0.012 

11.630±0.089 

10.236±0.154 

9.255±0.062 

9.180±0.057 

4.68 (4.35 ROB) 

10.5±2.0* 

97 

1.07 

4.50 (4.83 ROB) 

24 

10.33769 

-0.08369 


17.941±0.114 

15.592+0.027 

13.362±0.090 

12.502±0.116 

11.014±0.144 

9.411±0.217 

5.06 

0.9±0.3* 


1.19 

4.35 

25 

10.36273 

-0.05901 


17.617±0.085 

16.249±0.049 

14.468±0.177 

12.621±0.129 

11.510±0.133 

10.473±0.104 

6.5±0.5* 

<0.3* 


0.87 

3.97 

26 

10.23508 

-0.12661 




13.425±0.115 

11.594±0.128 

10.272±0.087 

9.854±0.104 

6.50±0.20* 

0.4±0.2* 



3.35 

27 

10.25199 

-0.07062 

18.642+0.089 

18.236±0.149 

16.330+0.053 

11.443 ±0.039 

9.736±0.045 

8.199±0.032 

7.340±0.038 

4.00±0.10* 

<1.0* 


1.86 

3.34 


Notes. The 24 iim magnitudes are from the MIPSGAL catalogue or from Robitaille et al. (2008; magnitude in brackets followed by ROB), or have 
been measured by us. A * indicates our own measurements. The luminosity of the sources (associated with GOlO.32-00.15) is estimated from 
their 70 fim flux, using the correlation established by Duham et al. (2008) and Ragan et al. (2012).The positions of the 6.7 Ghz methanol masers 
are from Green et al. (2010; positional accuracy ^CK.'d), the positions of the EGOs from Cyganowski et al. (2009; positional accuracy ~0'.'l). 
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Appendix A: More about G010.32-00.15 

A.1. Distance of G010.32-00.15 

The distance of this region is very uncertain, and different indi¬ 
cators (and authors) give distances covering the range 2 kpc to 
19 kpc. 

• Radio recombination lines give the velocity of the 
ionized gas: VLSR(H109a)=9.7 km s“^ (Reifenstein et 
al. 1970), V(H110a)=12.0 km s”! (Downes et al. 1980), 
V(H110Qf)=9.4 km s"^ in the central region (Araya et al. 2007). 
The molecular condensations present at the waist of the nebula 
have similar velocities in the range 11.3 - 13.5 km s“^ (NH 3 
emission lines; Wienen et al. 2012; see Table 4); see also the 
CO maps in Beuther et al. (2011; their figures 7 and 8 ). These 
velocities lead to a kinematic near distance of ~2 kpc or a far 
distance ~15 kpc. 

• H 2 CO absorption lines at 6 cm are often used to resolve 
the distance ambiguity. For distances within the Sun’s orbit, the 
radial velocity increases with distance from the Sun in the first 
quadrant, reaches a maximum at the tangent point and then de¬ 
creases to zero at the Galactocentric radius of the Sun. If the 
thermal continuum of the H ii region is absorbed by the H 2 CO 
molecules, the H 2 CO-absorbing cloud must be in front of the H ii 
region. H 2 CO absorption lines in the direction of GlO.31-0.15 
are observed at velocities of 11.1 km s“\ 22.3 km s“\ and 
38.4 km s“\ with faint components at 3.9 km s“\ 22.3 km s“\ 
and 47.3 km s-^ (Sewilo et al. 2004; HPBW 2.56'). Similar ve¬ 
locities were measured for the H 2 CO absorption lines by Wilson 
(1974) and by Downes et al. (1980). That absorption lines are 
observed at velocities higher than for the H ii region points to the 
far kinematic distance. A far distance was adopted by Georgelin 
& Georgelin (1976; 18.7 kpc), Sewilo et al. (2004; 15.0 kpc), 
and Du et al. (2011; 15.12 kpc). 

• The velocity at the tangent point is -150 km s“\ as indi¬ 
cated by the Hi emission observed in this direction (Kalberla 
et al. 1982), but the H 2 CO lines show no absorption near 
150 km s“^ (Sewilo et al. 2004). If GOIO.32-00.15 was at the far 
distance, we should observe absorption up to the maximum ve¬ 
locity in this direction. The H 2 CO absorption only extends up to 
47 km s"\ which could be due to a lack of H 2 CO material near 
the tangent point. But the same situation is observed with the 
atomic hydrogen: the Hi absorption extends only to 45 km s“^ 
in the direction of GOIO.3-00.1 (and of the two nearby Hii 
regions GOlO.2-00.3 and GOlO.6-00.4), whereas its emission 
extends to 152 km s“^ (Kalberla et al. 1982). Thus the Hii re¬ 
gion GOlO.32-00.15 cannot be at the far distance. This situation 
has been discussed by Wilson (1974) and Downes et al. (1980). 
They adopt a distance of -5-6 kpc, corresponding to the high¬ 
est H 2 CO absorption velocity. A distance of 6 kpc also has been 
adopted by Beuther et al. (2011) in their study of W31. 

• The nearby region GOlO.2-00.3 has velocities similar 
to our region, both for the ionized gas and for the associated 
molecular material (Wilson 1974; Downes et al. 1980; Araya et 
al. 2007; Beuther et al. 2011; Wienen et al. 2012), so these two 
regions probably lie nearby. Near-IR photometry and K-band 
spectroscopy of the exciting stars of GOlO.2-00.3 by Blum et 
al. (2001) indicate a distance of 3.4+0.3 kpc. 

We identified the exciting cluster of the Hii region in 
Sect. 6.1 and its main 05V-06V exciting star. The spectropho- 
tometric distance of this star is estimated to be 1.75 kpc, in good 
agreement with the near kinematic distance. And we have shown 
that this distance could not be greater than 2.5 kpc. 


How can we reconcile this maximum distance of 2.5 kpc for 
GOlO.32-00.15 with the H 2 CO and Hi absorption spectra? We 
note that this situation is not unique to GOlO.32-00.15. A paral¬ 
lax of a maser spot associated with GOIO.6-00.4 was measured 
by Sanna et al. (2014), and they find a distance of 4 . 95^^43 pc. 
Quireza et al. (2006) find a radio recombination line velocity of 
-1.1 km s“^ for GOIO.6-00.4, which leads to a kinematic dis¬ 
tance of -18 kpc (and no near kinematic distance). Similarly, the 
HII region G009.62-r00.19 has a parallax distance of 5. 15^^33 pc 
(Sanna et al., 2009) and near and far kinematic distances of 
-1 and -16kpc. Hofner et al. (1994) found Hi absorption up 
to 55 km s“^ for G009.62-1-00.19. All three sources are near 
a Galactic longitude of 10°, have LSR velocities <10 km s“\ 
have H 2 CO and HI absorption up to -50 km s“\ and have large 
discrepancies between their distances estimated using different 
indicators. The H 2 CO and Hi clouds that are seen in absorption 
up to -50 km s“^ must lie in front of the three sources. This 
implies that the velocity field cannot be described by approxi¬ 
mately circular orbits about the Galactic centre. The deviations 
from circular velocity here have been noticed by many previous 
authors and are probably the effects of the Galactic bar(s) (e.g., 
Rodriguez-Fernandez & Combes, 2008). Clearly, kinematic dis¬ 
tances in this part of the Galaxy are not accurate. 

A.2. Other candidate Class I YSOs in the vicinity of 
G010.32-00.15 

We discuss here the candidate Class I YSOs present in the vicin¬ 
ity of GOl0.32-00.15, which are not associated with it or whose 
association with the central H ii region is uncertain. 

• YSOs in the direction of C7: this clump does not belong 
to the complex; its distance is unknown, preventing us from es¬ 
timating its mass and density. 

The northern extension of C7 contains a small H ii region de¬ 
tected on the MAGPIS radio map and the faint candidate Class I 
YSO #14 (Fig. A.l (a)). YSO #14 has no detectable 24 yum and 
70 yum counterparts, and its nature is uncertain (a low-luminosity 
Class II or fiat-spectrum source). The flux of the UC Hii region 
is given in Table A.2; its exciting star has not been identified. Its 
size is -3'.'7 (beam deconvolved). The 70 yum source associated 
with this UC Hii region has a flux of 130+10 Jy. 

C7 contains a compact Hii region (see Table A.2) and the 
candidate Class I YSO #15 in a nearby direction. YSO #15 is 
bright at 24 yum ([8.0]-[24] in the range 4.30 mag to 4.47 mag) 
and at 70 yum, but at this last wavelength it is difficult to separate 
it from the extended emission linked to the compact H ii region. 
At 70 yum we measured a flux of 113+30 Jy for this Hii region 
plus YSO #15. A class II methanol maser was detected by Green 
et al. (2010) at less than 1" from YSO #15. The velocity of the 
maser, 39.0 km s“^ indicates that it is associated with C7 and 
not with the GOlO.32-00.15 bipolar nebula. Thus YSO #15 is 
clearly a Class I YSO embedded in Cl. 

A small Spitzer and Herschel-FACS diffuse region lies about 
25" west of YSO #15, source F in Fig. A.l. It has no detectable 
radio-continuum counterpart on the MAGPIS map (flux at 
20 cm <5 mJy). 

• Sources at the periphery of the bipolar nebula: they are 
identified in Fig. A.2. 

YSOs #16 and #17 lie at the extremity of the bottom lobe. 
YSO #16 presents a 24 pm and a faint 70 pm counterpart 
([8]-[24]=3.15 mag, hence an uncertain nature. Class II or fiat- 
spectrum YSO). YSO #17 has no detectable 24 pm and 70 pm 
counterparts so it is probably a low-luminosity Class II YSO. 
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Fig.A.l. Cl clump, (a) Colour composite image with the 24 yum, the 
8.0 jim, and the 3.6 jim emissions, respectively in red, green, and blue 
(linear units). Two candidate Class I YSOs, #14 and #15, are identified; 
the brightness of YSO #14 has been enhanced. The three blue arrows 
point to candidate AGB stars. The red contours are for the column den¬ 
sity (same levels as in Fig. 21 (b)). (b) The column density contours 
are superimposed on the MAGPIS radio map at 20-cm. The UC and 
compact H ii regions are identified, which have 24 fim and 70 fim coun¬ 
terparts. 


They are not associated with a condensation (the column density 
in their direction is ~2.5xl0^^ cm“^). YSOs #18 and #19 lie east 
of the southern lobe in a region where the column density is in 
the range 3 to 4x10^^ cm“^). They have a 24 jum counterpart but 
no detectable 70 yum one ([8]-[24] is in the range 2.53 mag to 
2.83 mag for YSO #18 and -3.01 mag for YSO #19). They are 
probably low-luminosity Class II YSOs. 

YSOs #24 and #25 lie in the direction of the northern lobe, 
in directions of intermediate column density (in the range 4 to 
5x10^^ cm“^). They have 24 yum counterparts and colours typical 
of Class I YSOs ([8]-[24]=4.35 mag for #24 and 3.97 mag for 
#25). YSO #24 has a faint 70 yum counterpart, and YSO #25 has 
none. 


YSO #27 lies in the direction of a nearby faint Hii region. 
We argue that it is probably associated with this H ii region and 
not with the bipolar nebula. (There is no velocity measurement 
to indicate the distance of this nearby nebula.) The association is 
uncertain for YSO #26, which lies in the direction of a filament 
west of the bipolar nebula (column density 4.8x10^^ cm“^). The 
colour of these sources indicates that they are probably flat- 
spectrum YSOs ([8]-[24]~3.35 mag for both). 
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Fig.A.2. Identification of the candidates Class I YSOs in the vicinity of 
GOlO.32-00.15. The underlying grey image is the Spitzer 8.0 pm im¬ 
age. The red contours are the column density contours (levels of 3,5,10, 
and 15 xlO^^ cm“^. The intensity of some YSOs has been enhanced.) 


A.3. Tables 

Appendix B: Apertures used to estimate the mass 
of the dumps 



Fig.B.l. Cl, C2, and C3 clumps in the field of G319.88 -l 00.79. The 
underlying grey image is the column density map. The apertures used 
to integrate the column density and estimate the masses are indicated 
(red contours; these apertures follow the level at half peak-values, after 
correction for a background emission). The green contours correspond 
to column densities of 1, 2, 4, and 6x10^^ cm"^ (as in Fig. 7 (d)). 
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Fig.B.2. Cl to C7 clumps in the field of GOlO.32-00.14. The under¬ 
lying grey image is the column density map (logarithmic units). The 
apertures used to integrate the column density and estimate the masses 
are indicated (green contours; the polygonial apertures follow the level 
at half peak-values, after correction for the background emission, the 
same for all clumps). The red contours correspond to column densities 
of 0.3, 0.5, 0.75, 1.0, and 1.5x1023 cm'^ (as in Fig. 15 (d)). 

Appendix C: Infrared dark clouds in the direction of 
the bipolar HII regions 



Fig. C.l. IRDCs detected by Peretto & Fuller (2009) in the vicinity of 
G319.88-^00.79. The underlying grey image is the Spitzer 8.0 pm im¬ 
age. The red contours correspond to column densities of 1, 2, 4, and 
6x10^2 cm-2. 



Fig. C.2. IRDCs detected by Peretto & Fuller (2009) in the vicinity of 
GOlO.32-00.15. The underlying grey image is the Spitzer 8.0 pm im¬ 
age. The red contours correspond to column densities of 0.5, 1.0, 1.5, 
and 2.0x10^2 cm"^. 
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Table A.l. Exciting cluster of GOlO.32-00.15. Data extracted from the UKIDSS catalogue 
except for star #1. 


N 

1 

b 

j 

H 

K 



n 

O 

(mag) 

(mag) 

(mag) 

(") 

1^ 

10.31798 

-0.15185 

11.631+0.043 

9.694+0.055 

8.613+0.039 

0.00 

1' 



11.9 


9.0 

0.00 

2 

10.31773 

-0.15134 

13.662+0.002 

11.915+0.001 

10.875+0.001 

2.16 

3 

10.31872 

-0.15123 

14.065+0.002 

12.180+0.001 

11.034+0.001 

3.47 

4 

10.31812 

-0.15058 

17.175+0.024 

14.874+0.007 

13.578+0.005 

4.79 

5 

10.31686 

-0.15266 

17.247+0.026 

14.783+0.007 

13.405+0.004 

4.99 

6 

10.31945 

-0.15219 

14.286+0.002 

12.392+0.001 

11.316+0.001 

5.61 

7 

10.31958 

-0.15174 

15.372+0.005 

13.337+0.002 

12.005+0.001 

5.83 

8 

10.31849 

-0.15348 

14.196+0.002 

12.079+0.001 

10.789+0.001 

6.53 

9 

10.31614 

-0.15226 

17.224+0.025 

15.190+0.009 

14.049+0.007 

6.81 

10 

10.31641 

-0.15316 

15.404+0.005 

13.481+0.002 

12.392+0.002 

7.39 

11 

10.32038 

-0.15154 

16.169+0.010 

14.141+0.004 

12.969+0.003 

8.78 

12 

10.31986 

-0.15340 

17.306+0.027 

14.910+0.007 

13.472+0.004 

9.28 

13 

10.31880 

-0.14947 

12.105+0.001 

11.535+0.001 

10.750+0.001 

9.32 

14 

10.32000 

-0.15004 

15.917+0.008 

15.075+0.008 

14.595+0.011 

9.81 

15 

10.32061 

-0.15105 

16.503+0.013 

14.435+0.005 

13.190+0.003 

9.92 


^ d is the projected distance from the exciting star. 
^ 2MASS photometric data. 

^ photometry from Bik et al.(2005). 


Table A.2. Ultracompact Hii regions near GOlO.32-00.15. 


Name 

Radio fiux, A 
(mJy, cm) 

Diameter 

(") 

Reference 

C3 




GOlO.301-00.147 

654.6, 6 


Becker et al. (1994) 

GOlO.30058-00.14734 

426.2, 20 


Helfand et al. (2006) 

GOlO.3009-00.1477 

631.4,5 

5.2 

Urquhart et al. (2013) 

GOlO.321-00.147 

320.7,20 

4.9X4.2 

White et al. (2005) 

GOlO.321-00.147 

628.0,6 

5.5X4.1 

White et al. (2005) 

Cl extension 




GOlO.321-00.233 

67.9, 6 


Becker et al. (1994) 

GOlO.32029-00.23323 

56.7, 20 


Helfand et al. (2006) 



3.7* 

MAGPIS 

Cl extension 




Source B 

8% 20 

unresolved* 

MAGPIS 


Notes. Our own measurements are indicated by an *. The radio continuum source G10.321-0.147 is a UC Hii region as defined by Kurtz (2005): 
its size is ~0.045 pc<0.1 pc, its emission measure is ~10^ cm“^ pc, its density is 1.9x10^ cm >10"^ cm“^, its mass is ~2xl0“^ M©. The nature 
of the two other radio sources is uncertain because their distance is unknown, but they have small sizes. 
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